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ABSTRACT 

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the 
galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing 
Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected elliptic¬ 
ity components of galaxy matter haloes and their embedded galaxies, /h, split by galaxy type. 
We then compare our observations to measurements taken from the Millennium Simulation, 
assuming different models of galaxy-halo misalignment. Using the Millennium Simulation 
we verify that the statistical estimator used removes contamination from cosmic shear. We 
also detect an additional signal in the simulation, which we interpret as the impact of intrinsic 
shape-shear alignments between the lenses and their large-scale structure environment. These 
alignments are likely to have caused some of the previous observational constraints on fy to 
be biased high. From CFHTLenS we find /h = —0.04 ± 0.25 for early-type galaxies, which 
is consistent with current models for the galaxy-halo misalignment predicting fy 0.20. For 
late-type galaxies we measure /h = 0.6!) JJ'jy from CFHTLenS. This can be compared to the 
simulated results which yield fy ~ 0.02 for misaligned late-type models. 

Key words: gravitational lensing: weak - galaxies: haloes. 


1 INTRODUCTION 

In the standard cosmological paradigm, galaxies, groups, and 
clusters are embedded in large haloes of - mostly dark - 

* E-mail: schrabba@astro.uni-bonn.de 


matter. Numerical simulations of cosmic structur e formation 
predi ct that these haloes are roughly triaxial (e.g. Ijing & Sutol 
l2002ll . and that their average density profile s close ly follow the 
Nava rro-Frenk-White profile (NFW, iNavarro. Frenk & White! 
1199(i 11997b . In projection, these should approximately appear 
elliptical. This prediction can be tested observationally on the 
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scales of galaxies and clusters. One approach is to constrain halo 
shape s via the use of b aryonic tracers s u ch as satellite galaxies 


(e.g. iHolmberd 19691: iBrainer 3 20051: 


Yang et af||2006 : Azzaro etall 12(107 


Libeskind et al. 


g alaxies 
1. 2005: 


MetZjJCrou^a&Jerjenl 


Nierenberg et al 


20071: iFaltenbacher et al^ 

Okumura, Jing & Li 12009 


200 71 120091: i Bailin etal 

lAgustsson & Brainer j 


200S; 


2010 


__ ,201 ill, the distribution of stella r velocity 

l Oiling & Merrifielcf 2000h. satellite tidal st r eams dlbata et all 
l200ll: iLuxetalJ 1201A IVera-Ciro & Helmil 120131). HI ga s 


1 Baneriee & Jog TioOSl : (yBrier^Freeman&v an der Kruitl l2010h . 

or planetary nebulae dNapolitano et al. 201 lh . Such constraints 


can be compared to the output of hydrodynamical simulations 
that aim at modelling galaxy formation, linking the baryonic and 
dark matter components, and which can provide predictions on the 
relative (mis-) alignment of galaxies, their dark matter hosts, and 


the surrounding large-scale structure, as well 

as on the impact of 

the baryons on halo shapes (e.g. Bailin et al. 

l2005i; Knebe et alj 

20id;IVera-Ciro et al. 

2011; Tenneti et al. 2014 

; Laigle et al. 20151: 

Debattista et al. 20151: 

Velliscig et al. 2015|). 


As an observational alternative, gravitational lensing probes 
the total projected mass distribution directly without relying on vis¬ 
ible tracers. Strong gravitational lensing can provide information 
on t he inner shapes of the mass distribut ion in massive clusters 
(e.g. iLim ousin et~al] 12013 b and gala xies ( Ivan de Yen et al. M; 


iDutton et al, 2011 : Suvu et al]|2012h . but at these scales baryons 


have a non-negligible influence. At larger scales, constraints on the 
projected mass distribution can be obtained with weak gravitational 
lensing, which probes the coherent distortions imprinted onto the 
observed shapes of background galaxies by the tidal grav itational 
field (e.g. lBartelmann & Schneidej200ll : ISchneideil2006l) . Even at 
these scales baryons contribute to the total mass distribution probed 
by weak lensing (e.g. in the form of satellite galaxies), but the dark 
matter is expected to dominate. 

Weak gravitational lensing has been successfully used in a 
wide range of applications, including for example cosmic shear 


(e.g. 

Schrabback et al. 

2013; Kilbinger et al. 2013; Hevmans et al. 

2013 

: Huff et al. 2014; 

Kitching et al. 20141: Simon et al. 2015f the 

mass calibration of galaxy clusters (e.g. von der Linden et al. 20141; 

Kettula et al. 2015]; Ford et al. 2015: Hoekstra et al. 20150, and 

constraints on the azii 
matter haloes (e.g. 

nuthally-averaged mass profiles of galaxy 

Vlandelbaum et al. 2006c; van Uitert et al. 

2011 

; Leauthaud et al. 2012j; Velander et al. 2014; Hudson et al. 


2015,1 


In principle, weak lensing is also sensitive to halo elliptic- 
ity as the gravitational shear at a given radius is larger along 
the direction of the major axis of the projected halo compared 
to the projected minor axis (e.g. iNataraian & Refregieil l2000l : 
IBrainerd & Wrightl 120001) . For very massive clusters, the weak- 
lensing signal is strong enough to provide individual halo ellip - 
ticity constraints ICorless. King & Clowell2009l : l0guri et al.l20ld) . 
Constraints were also obtained in stacked analyses of larger clus¬ 
ter samples fevans & Bridlel2009l:l0guri et allioll) . For less mas¬ 
sive, galaxy-scale haloes, the signal can still be detected statistically 
by stacking very large samples. 

A measurement of this effect through stacking requires one 
to align the shear field around all lenses prior to stacking, such 
that the major axes of all projected haloes are aligned. The ori¬ 
entations of the (mostly dark) matter haloes, however, are not di¬ 
rectly observable. One approach which has been used instead for 
weak lensing halo shape studies, is to approximate the orienta¬ 
tions of the projected haloes by the orientations of their galaxy im- 


ages 1 Hpekxtra^jfee&Cjladdgrs 2004 : Mandelbaum et al] l2006al : 


IParker et ah 20071 : van Uitert etalT 20121 . Also, it has typically 

been assumed that, on average and in projection, more ellipti¬ 
cal lenses are hosted by more elliptical halo es, as supp orted by 
simulations for early type galaxies dNovak et al J 12006 ). A key 
parameter which is then extracted is the average aligned ellip- 
ticity ratio between the ellipticities of the projected halo and 
the observed lens light distribution /h = (cos(2A^h, g )|eh|/|e g |), 
where the averaging typically includes a weighting scheme that 
depends on |e g |. Here, A0h, g is the angle between the ma¬ 
jor axis of the projected galaxy light distribution and the ma¬ 
jor axis of its projected matter halo. Thus, only in the case 
of perfect alignment (A^h.g = 0), fh reduces to the actual el- 
lipticity ratio. However, in practice one expects a considerable 
random misalignment between the observed shapes of galaxies 
and matter haloes, as su ggested both by numerical simulations 
of galaxy form a tion (e.g. OkamotoetalJ 2005); CrainrtaL 20091 
Be^et_alJ|20d(i yDeasom et al .11201 ll : iDong et al] 20141 : IWang et all 


2Q14; iTenneti et al.ll2014h and observations t hat approximate mat¬ 

ter ha loes via the distribution of satellites (e.g. lOkumura. Jing & Lil 
l2009ll . This should substantially reduce /h and wash out the halo 
shape signature, making it difficult to detect observationally ( Betti 
l2012ll . It is important to test this prediction observationally, both to 
improve our understanding of galaxy formation, but also to inform 
models of intrinsic galaxy alignments. Such alignments of galaxies 
with their surrounding mass distribution are an important physi cal 
contaminant for cosm ic shear studies (e.g. iHevmans et alj|2013l : 
Ijoachimi et al.ll2013bl). Of particular con cern are shape-shear in¬ 
trinsic alignments 1 Hirata & SeliaklljQoil) . Here, the ellipticities of 
foreground galaxies are aligned with their surrounding large-scale 
structure, which lenses the background sources. So far, most con¬ 
straints on shape-shear alignments come from studies investigating 
the align ment of galaxies with the i r surrounding galaxy distribu¬ 
tion (e^ MandelbaumetaL 201T; Joachimietal 201 l]^ LJ et al J 


120131 : IZhang et alj 20131 : Singh, Mandelbaum & Moref 20151) . In- 
terestingly, the aligned halo shape signature we want to extract di¬ 
rectly contribut es to the shape-shea r intrinsic alignment signal at 
small scales jBridle & Abdallall2007 ;t. 

In addition to the expected small signal, there are further ob¬ 
servational challenges: While a potential additional alignment of 
lenses and sources does not affect the azimuthally-averaged galaxy- 
galaxy lensing signal, this is no longer the case for the most simple 
estimator of the anisotropic halo ellipticity signal. Such a source- 
lens alignment can, for example, be introduced by incomplete re¬ 
moval of instrumental signatures such as the point-spread function 
(PSF). In addition, gravitational lensing by structures in front of 
the lens causes extra alignment. The latter is commonl y ref erred 
to as cosmic shear, or sometimes as multiple d eflections (IBrainerd 
120101) . For example. fHowell & Brainerd J20 1 (j ) study the impact of 
cosmic shear on halo ellipticity constraints and conclude that obser¬ 
vational estimates of the ratio of the shears along the lens major and 
minor axes would need to be compared to Monte Carlo simulations 
for interpretation. Ho wever, they do not consider the modif ied es - 
timator introduced by Mandelba um et al] d2006aL hereafter lM0(il) . 
which leads to a cancellation of such spurious signal at the relevant 
(small) scales, and which we employ in the current study. 

Previous observational constraints on halo ellipticity from 
weak lensin g are_somewhat inconclusive: ^)n the one hand, 
iHoekstra. Yee & Gladdersl d2004lf and IParker et aP (120071 1 find in¬ 
dications for positive /h using magnitude-selected lens samples 
and simple esti mator s th at do not correct for sy stematic shear. On 
the other hand. iMO^ and lvan Uitert et alj ( 120120 separate lenses by 


© 2015 RAS, MNRAS POQ.fTHTTl 

















































































































































































































CFHTLenS: Halo ellipticity constraints 3 


colour and correct for systematic shear , but do not detect signifi¬ 
cantly non-zero fh. For example. lM06l find fh = 0.60 ± 0.38 for 
red and fh = — lA^'l for blue lenses using d ata fr o m the Sloan 
Digital Sky Survey (SDSS). Very recently. I Clampitt & Jainl l l2015l) 
reported a significant detection of the signature of halo ellipticity 
around luminous red galaxies, also employing SDSS data. 

In this study we use weak lensing data from the 
Canada-France-Hawa i i Telescope Lensing S u rvey (CFHTLenS , 
He^mans_etaL 20121 : iHildebrandt et alj 1 20121 : lErben et all 1 201 3l : 


Miller et al.1120 1 3h ~to derive updated constraints on the ellipticity 


of galaxy matter haloes from weak lensing. These data allow us to 
constrain the signal for galaxies subdivided into bins of photomet¬ 
ric type and stellar mass. The former division allows us to approx¬ 
imately separate the lens sample into early and late type galaxies, 
which have different predictions for the expected weak lensing halo 
shape signal. The latter division optimises the total measurement 
signal-to-noise as stellar mass acts as a proxy for halo mass and 
therefore the signal strength. 

In addition, we study a simulated weak lensing halo shap e sig¬ 
nal based on the Millennium Simulatio n d.Snrinsjel et alj2005ll . em¬ 
ploying the ray-tra cing analysis fro m Hilbert et alj J2009h and lens 
shapes computed in ljoachimi et alj d201 3al lblf. On the one hand, this 
provides a sig nal predicti on given current galaxy-halo (mis-) align¬ 
ment models d Bettll20 ill that we can compare the CFHTLenS re¬ 
sults to. On the other hand, it allows us to study the impact of cos¬ 
mic shear on the measurement and check for deviations from the 
simple model prediction of isolated elliptical NFW haloes. 

This paper is organised as follows: Sect.[2] describes the for¬ 
malism of the halo ellipticity measurements using weak lensing, 
and the verification tests which we have conducted using a simple 
simulation. In Sect.[3]we summarise properties of the CFHTLenS 
data, discuss the selection of lens and source galaxies, and present 
the measured shear signal and constraints on the aligned ellipticity 
ratio fh- We present our analysis of the simulated data based on the 
Millennium Simulation in Sect.[4] We then discuss our results and 
conclude in Sect.[5] In addition, we present a consistency check for 
shape measurements of background sources in the vicinity of bright 
foreground lenses using image simulations in Appendix lAl 

For the computation of angular diameter distances we by de¬ 
fault assume a flat ACDM cosmology with D m = 0.3, Da = 0.7, 
Hq = 70/i7o km/s/Mpc, and hjo = 1, which is consistent with 


the best-fittin g cosm ological parame ters f rom both WMAP9 
dHinshaw et all 120131 ) and Planck l lPlanck Collaboration et alj 
|2015i) at the ~ 1 — 2cr level. The only exception is our analy¬ 
sis of the simulated data from the Millennium Simulation, for 
which we use the input cosmological parameters of the simulation 
D m = 0.25, Da = 0.75, Ho = 73/173 km/s/Mpc, and (173 = 1 
iSpringel et al.l20C)3) . All magnitudes are in the AB system. Stellar 
masses M, are given in units of solar masses Mg. 


2 METHOD 


2.1.1 Constraining the isotropic galaxy-galaxy lensing signal 

In weak lensing studies the shape of a galaxy is typically described 
by the complex ellipticity 

e = ei + ie 2 = |e|e 2l<A . ( 1 ) 

For the ellipticity definition employed here (see lMiller et al J201 3l Pl 
and the case of an idealised source with elliptical isophotes, the 
absolute value of the ellipticity is given by |e| — (a — b)/(a + b). 
Here, a and b are the major and minor axes of the ellipse, while 4> 
corresponds to the position angle of the major axis from the x-axis 
of the coordinate system. In this definition the ellipticity of each 
background galaxy provides an unbiased but very noisy estimate of 
the reduced gravitational shear g in the direction of the source 

E(e) = g = ~ 7 , (2) 

1 — « 


where E indicates the expectation value, while 7 and « denote the 
shear and convergence. In principle, all structures along the line of 
sight from the source to the observer c ontribute to the net shear and 
convergence (see e.g. lSchneideifeOOfii) . However, in galaxy-galaxy 
weak lensing analyses, such as the study presented here, one corre¬ 
lates the shear inferred from background sources with the positions 
of foreground lenses. In this case, only structures at the lens redshift 
contribute to the net signal (other structures add noisfl), allowing 
us to express the correlated convergence k = E/E c as the ratio 
between the surface mass density E and the critical surface mass 
density 


E c 


_c^_1_ 

47tG D\/3 


(3) 


Here, c denotes the speed of light in vacuum and G the gravitational 
constant. The geometric lensing efficiency f3 is defined as 


P 


max 



(4) 


D s , D \, and D\ B indicate the physical angular diameter distances to 
the source, to the lens, and betw een lens and source, respectively 
(note that the equations in [M06i are expressed for comoving dis¬ 
tances instead). 

Galaxy-galaxy weak lensing analyses study the stacked shear 
field around foreground lens galaxies. For this it is useful to decom¬ 
pose the shear and correspondingly the ellipticities of background 
galaxies into the tangential component and the 45 degrees-rotated 
cross component 


et = — ei cos 2d — e 2 sin 20 , (5) 

e x = +ei sin 26 — e 2 cos 26 , ( 6 ) 


where 6 is the azimuthal angle with respect to the lens position as 
measured from the x-axis of the coordinate system used. 

The majority of the previous galaxy-galaxy weak lensing anal¬ 
yses have only studied the profile of the azimuthally-averaged tan¬ 
gential shear 7 t, which relates to the differential surface mass den¬ 
sity profile AE(r) = E(< r) — E(r), where E(< r) is the mean 


2.1 Formalism 

The methodology of our analysis large ly follows the approach and 
notation that was intr oduced by M06 and additionally applied in 
Ivan Uitert et al.1 d 2012 h . It allows for the correction of spurious sig¬ 
nal originating from cosmic shear or instrumental distortions. 


1 This ellipticity definition is often r eferred to as e in the literature. Here, 
we denote it as e to be consistent with lMiller et al.l JdO 1 3| ). 

2 For a weak lensing halo shape analysis foreground structures in front of 
the lens cause an extra alignment of sources and lenses, which does in¬ 
troduce spurious signal for simple estimators. However, this signal is ac¬ 
counted for via the formalism explained further below in this section. 
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convergence within radius r, as AE = 7 tE c . We estimate the dif¬ 
ferential surface density from the source galaxy ellipticities as 


AE(r) = 


Ei WiE Cii (et,jE Cj j) Ei W *E,i e t,i 


E t ^EE 2 




(7) 


where we sum over lens-source pairs in an annulus around r. In this 
section we indicate estimators with a hat for clarity, but we drop it 
in the subsequent sections when presenting results. In our analysis 
we employ inverse-variance weights Wi for the source shape esti¬ 
mates. These weights account for both measurem ent noise and the 
intrinsic ellipticity distribution jMi 11 er et al J 20 1 3l ). We use uniform 
source shape weights for simulated data. The main reason for con¬ 
ducting the analysis in terms of AE, which is a re-scaled version 
of the shear, instead of the shear directly, is to adequately account 
for the redshift dependence of the weak lensing signal. 

We note that the tangential ellipticity components et of 
sources provide estimates for the tangential component of the re¬ 
duced shear gt, while AE is defined in terms of the tangential 
component of shear 7 t. In galaxy-galaxy weak lensing, typically 
| 7 | 1 and |k| 1. Hence, many studies have typically approxi¬ 

mated the reduced shear with the shear as indicated in Here we 
implicitly account for the difference when fitting the azimuthally- 
averaged tangential shear profiles, as we find a small but non- 
negligible impact for our most massive lenses. When studying the 
anisotropy in the shear field, as detailed below, we however ignore 
reduced shear corrections as they cancel out to leading order. 

In our analysis we fit the isotropic part of the measured 
shear profile with a n NFW shear profile prediction according to 
Iwright & Braineixl feood t in order to constrain r 2 ooc, the radius 
corresponding to a mean over-density that is 200 times the critical 
density at the lens redshift, from the data itself. For this , we employ 
the m ass-concentration relation of NFW haloes from iDuffv et all 

d 2008 l) . 


2.1.2 Constraining the anisotropic galaxy-galaxy lensing signal 

The formalism to study the anisotropic w eak lensing shear field 
around elliptical lenses was introduced by INataraian & Refregien 
( I 2 OOOI) for the case of an elliptical isothermal sphere, and f urther 
developed and generalised for other density profiles in lVIOfil . Here 
we largely follow the notation from lM06l . and introduce a few ad¬ 
ditional quantities;_ 

Similarly to lM06l . we model the stacked and scaled tangential 
shear field as a combination of an isotropic profile AEj SO (r) and 
some azimuthal variation as 

AE mode i(r, Ad) = AEi SO (r) [1 + 4/ re i(r)|e h ,a| cos(2A6>)] . 

( 8 ) 

Here, Ad denotes the position angle with respect to the major axis 
of the lens galaxy. We do not know the orientations of the matter 
halo ellipticities eh on the sky. Thus, we have to approximate them 
with the orientations of their corresponding galaxy ellipticities e g 
when stacking the anisotropic shear field. Accordingly, our analysis 
is only sensitive to the average component 

|eh,a| = (cos(2A0h,g)|eh|) ~ (cos(2A(/> h , g ))|eh| (9) 

of the halo ellipticity that is aligned with the galaxy ellipticity, 
where A^>h, g indicates the misalignment angle. Here we average 
over the misalignment distribution, which we assume does not de¬ 
pend on | eh | - Following lM06l . we make the assumption that the 
absolute value of the halo ellipticity is proportional to the absolute 


value of the galaxy ellipticity 

|e h |=/h|e g |. ( 10 ) 

While there will be deviations from this assumed linear scaling in 
reality, it pro v ides a reaso nable approximate weighting scheme (see 
also lvan Uitert et alkoil who explore additional schemes). 

In equation J 8 ). / re i(r) describes the relative asymmetry in the 
shear field for an elliptical halo of ellipticity |eh, a |. It depends on 
the assumed density profile and n eeds to be computed numerically 
for non-power-law profiles (see lMOblPl . To recover the notation of 
[mqE we define 

fh = /h(cos(2A0 h , g )} = -p4l(cos(2A0h, g )), (11) 

\ e B\ 

f{r ) = /relM/h- (12) 

Then, equation ([Sji reduces to 


AE mode i(r, AO) = AE iso (r) [1 + 4/(r)|e g | cos(2A0)] . (13) 

lM06i show that the joint solution for the estimators of the isotropic 
and anisotropic shear field components is given by 


AEis 0 (r) 


/(r)AEiso(r) 


E i et :i 

Ei Wi E,i 

E i tw»E“j e t ,»|e g ,i| cos(2A0j) 
4 Ei w *E~ l 2 |e g , i | 2 cos 2 (2A6>;) ’ 


(14) 

(15) 


where the summation is again over lens-source pairs in a separation 
interval around r. Note that the factor 2 differen ce i n equations 
( 113b and ( 115b compared to equations (4) and ( 6 ) in M06l originates 
from the different ellipticity definition used by lM06l To ease the 
comparison to 1 m 06| we decided to not rescale /(r), but rather to 
write out the factor 2 difference explicitly. 

In practice, ([[5} is not a useful estimator for constraining halo 
ellipticity as it is susceptible to a systematic signal if the ellipticities 
of lenses and sources are aligned because of an additional effect. 
This could arise from instrumental systematics such as imperfectly 
corrected PSF anisotropy, but also from cosmic shear by structures 
in front of the lens. This can easily be understood: For example, an 
intrinsically round lens (with an isotropic halo) would appear ellip¬ 
tical because of this extra shear or systematic. Sources would also 
have an extra shear component parallel to the lens ellipticity. In the 
coordinates defined by the observed lens ellipticity this appears as 
an increased shear along the lens minor axis and a decreased shear 
along the lens major axis. Accordingly, this would be interpreted as 
an “anti-aligned" halo with /h < 0 . 

To cancel this systematic contribution. lM06l suggest to include 
an additional term in the estimator that is based on the ellipticity 
cross component e x ©, and which is given by 


/ 45 (r)AE iso (r) 


Ei WiEjex.ileg.il sin(2A0O 
4 Ei«'iE- l 2 |e g ,i| 2 sin 2 (2A6> i ) ' 


Equations ( 113b and d! 6 b obtain nearly equal contributions from sys¬ 
tematic effects aligning the lens and source ellipticities, as long as 
the shear correlation function £_(r) = { 7 t 7 t — 7 x 7 x)(f) is suffi¬ 
ciently small, where 7 indicates the additional “systematic” shear. 
This is the case at the relevant small r_(see our test for cosmic 
shear in Sect.[4] and the discussion in lM06l) . Hence, the estimator 


3 Our tests conducted in Sect. l2.2l indicate that there is only a weak depen¬ 
dence of /rel(r) on the halo ellipticity itself, which can be ignored for the 
expected halo ellipticities. 
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Uir) - Us(r)) AEi S o(r) can be used to probe halo ellipticity free 
from the systematic contribution. 

Importantly, / 45 (f) also contains signal from the flattened 
halo, so that both f(r) and /45(f) need to be modelled. Here we 
scale the model for J' 45 ( r ) in correspondence to equation d!2t as 

/45(f) = /rel,45(f)/h ■ (17) 

Im06I compute predictions for / re i(r) and / re i,45(f) for several el¬ 
liptical density profiles numerically. We restrict our analysis to el¬ 
liptical NFW profiles as they provide good fits to the isotropic shear 
signal of the CFHTLenS data within the cons idered radial range, 
see Sect. l3.5.TI For NFW density profiles lM06l compute model pre¬ 
dictions in terms of r/r s , where r s is the NFW scale radi us (w e 
linearly interpolate between the discrete values provided bv lMObl) . 
Accordingly, / re i(f) and / re i,45(f) can be computed for arbitrary 
NFW prof iles an d halo masses. We verify the numerical predic¬ 
tions from lM03 in Sect. 12.21 where we find very good agreement 
with lM06l for / re i(f). However, our analysis indicates an opposite 
sign for / re i,45(f) compared to the iMOq modelfl (see Sect. l2.2t . 
Hence, while we use the model interpolation scheme from IM061 
we implicitly use an opposite sign for / re i,45(f)- 


2.1.3 Estimating fh 


It is a primary goal of our analysis to constrain the aligned elliptic¬ 
ity ratio /h. From the estimators in equations G3- CL3. and m, 
and from the model descriptions in equations ( 1121) and ( 1171 . one 
could be tempted to define an estimator as 


with 

V(r ) 

x(r) 


/biased ( r ) 


yjr) 

X (f) 


_1_ 

/rel(f) - /rei,45(f) 

AEiso(f) . 


(/-/45)AE iso (f), 


(18) 


(19) 

( 20 ) 


However, the estimator ( 1181 would be biased due to the occurrence 
of the noisy x(r) in the denominator. In addition, we want to com¬ 
bine the estimates from the different radial bins, and possibly also 
different galaxy samples. One solution could be to divide by the 
best-fit model for AE; so (f), e.g. obtained from an NFW shear pro¬ 
filefitHowever, here we ado pt the alternat ive app roach suggeste d 
by iMOfil which is based on IBlissI d 1935allbh and iFielletl Il954f) : 
We are interested in the ratio m = y/x of two random variables. 
In our case, m corresponds to fh assuming / re i(f) and / re i,45(f) 
model the relative asymmetry in the shear field well. We assume 
that y and x have a Gaussian distribution, which is a reasonable 
approximation in galaxy-galaxy lensing given the dominant shape 
noise. We have multiple estimators yi, x t from the different ra¬ 
dial bins. Also, in some cases we want to combine the constraints 
from multiple galaxy samples originating from different redshift or 
stellar mass bins, which provide additional yi, Xi. For each i, the 
quantity yi — rrtXi is a Gaussian random variable drawn from an 
N(y, = 0, a 2 = wf 1 ) distribution with wf 1 = a~. + ma~.. Ac¬ 
cordingly, the summation 


E iWiiVi-mxi) f 1 \ 
~ V’E^J 


( 21 ) 


4 We note that different sign definitions exist in the literature for equation 
{6j. which might have led to an inconsistent model prediction derived in 

lM0d 


over all measurements is also a Gaussian random variable, where 
the distribution is taken at fixed m. We then determine frequentist 
confidence intervals at the Za level as 

-Z T,jWi(yi -mxj) z 
y/EEE Eitfi v 

By stepping through a grid in m we identify the best-fitting value 
that provides the desired estimator 

fh = m(Z = 0), (23) 


as well as 68 per cent confidence limits m(Z = ±1). At the best¬ 
fitting m(Z = 0) we also compute a reduced \' 2 of the fit as 


X 2 /d.o.f. 


EU - rn(Z = 0)f i ) 2 

n — 1 


(24) 


This approach assumes that off-diagonal covariance elements can 
be neglected. We find that this is indeed the case (see Sections [3.5.1l 
and l4Tt . 


2.2 Shear-field of an elliptical NFW halo 

In order to test our analysis pipeline and the model predictions from 
iMOd we first analyse a simplistic simulation containing a single, 
elliptical NFW halo. A more detailed simulation is presented in 
Sect.|3 using data from the Millennium Simulation and including 
more realistic galaxy samples, galaxy misalignment, and cosmic 
shear. 

For the basic test prese nted here, we generate an isotropic 
NFW convergence n profile iWright & Brainerdl 120001) on a fine 
4096 2 grid (10' x 1 0') for a M 2 q 0c = 10 12 M e /h 70 halo at 
z\ = 0.3 assuming the lDuffv et al.l 12008 ) mass-concentration rela¬ 
tion and sources at z s = 1. We shear the k profile along the tc-axis 
of the grid with various strengths to make it elliptical (eh,i > 0). 
Then the major and minor axes a and b of k iso-density contours 
fulfil |eh| = (a — b )/(a + b) . Ne xt, we co mpute the corresponding 
shear field using the lKaiser & Squires! (1 1 9931) formalism. To test our 
fitting procedure with a regular covariance matrix we generate 50 
shape noise realisations with a very small ellipticity dispersion (ar¬ 
bitrarily chosen cr e = 0.004), which were added to the shear field. 
For this test we simply use the halo ellipticity as the lens ellipticity, 
hence /h = 1. 

We plot the measured isotropic and anisotropic shear profiles 
for an |eh| = 0.2 halo in Fig. |T| The anisotropic signal (open sym¬ 
bols in the_bottom panel) is well described by the best-fitting model 
using the 1 mQ 6| prescription once we account for the sign error for 
/45 (see Sect. I TU Here we accurately recover the input halo el¬ 
lipticity (f h = 1.002±°;°°|). For comparison, we also show the 
measured signal if we artificially add a constant (scale indepen¬ 
dent) spurious alignment of the lens and sources corresponding to 
a shear of 72 = 19(0, to roughly illustrate the effect of residual PSF 
anisotropy or cosmic shear (cross- and star-shaped symbols in the 
bottom panel of Fig.|T}- While the /AE and —/4s AE components 
are individually disturbed, their combination (/ — /4s)AE is un¬ 
affected by the systematic shear as expected. 

5 If we always orient the lens for this test along the x-axis (ei only), we 
need to apply a systematic shear along the uncorrelated 45 degrees-rotated 
field diagonal ( 72 ). We verified that we obtain identical results with a sys¬ 
tematic shear along a random direction if we randomly rotate the lens and 
its shear field prior to applying the systematic shear. This would resemble 
reality more closely with random intrinsic lens orientations with respect to 
the orientation of the PSF or the cosmic shear field. 
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Figure 1. Test of the models and our analysis pipeline using a simple 
NFW k profile for a M 20 OC = 10 12 Mq//i7o halo at z\ = 0.3 and z s = 1 
sources. The halo was sheared on a grid to an elliptic ity eh = 0.2 and 
the co rresponding shear field was computed according to lKaiser & Squired 
<1993l) . The top and bottom panels show the isotropic and anisotropic shear 
profiles, respectively. For the anisotropic signal (note the scaling by r) the 
open blue squares, red hexagons, and black circles show /AS, — / 4 s AE, 
and (/ — / 4 s) AS, r espectively. The curves show corresponding model 
predictions from lMO^ for eh = 0.2 and /h = 1, where we have swapped 
the sign for the —/4s model (see Sect. Elq. The additional symbols show 
the same measured quantities, but with an extra constant 1 % shear added 
to both lens and source ellipticities (45° rotated from the lens orientation). 
The vertical dotted lines indicate the fit range. 


CM 



Figure 2. Test of the h alo ell ipticity recovery using the linearised elliptical 
shear field model from lMO^ (corrected for the sign inconsistency for / 4 s). 
Here the same type of elliptical NFW models was used as in Fig.|T] but for 
a range of halo ellipticities eh- Perfect recovery corresponds to the dotted 
/h = 1 line. Deviations for typical halo ellipticities are only a few per cent. 


We repeat this analysis with various halo ellipticities eh and 
plot the recovered /h(eh) in Fig.[2| This shows that the linear rela¬ 
tion assumed in equations o and is a good approximation, 
but leads to a slight overestimation of /h for very elliptical haloes. 
However, as the deviations from /h = 1 are small for typical ex¬ 
pected halo ellipticities^ of |eh| ~ 0.2, and given the large statis¬ 
tical uncertainties of CFHTLenS, we ignore this deviation in our 
analysis. 


3 OBSERVATIONAL CONSTRAINTS FROM CFHTLENS 


3.1 The CFHTLenS data 

To constrain the aligned projected ellipticity ratio fh observation- 
ally, we employ weak lensing data from the Canada-France-Hawaii 
Telescope Lensing Survey (CFHTLenS). It has an effective area of 
154 deg 2 imaged in the ugriz broad band filters using MegaCam 
on CFHT with a 5cr limiting magnitude in the det ection z-band 
and 2" apertures of zab ~ 24.5 — 24.7 dErben et al.ll2013l) . These 
data were obtained as part of the Wide-component of the Canada- 
France-Hawaii Telescope (CFHT) Legacy Survey which completed 
observations in early 2009. 

The CFHTLenS team has re duced the imaging data using 
the THELI pipeline dErben et al.l 120091 . l2013h . measured PSF- 
corre ct ed ga la xy shapes in t h e z-ba nd using lensfit dMiller et alj 
I2007L 120131 : iKitching et al.l l2008h . and estimated photomet¬ 
ric redshifts (ph oto-zs) from t h e ugriz d ata e mploying the 
BPZ algorithm |Benftez| I 2 OO 0 I : ICoe et all 120061) as detailed 
1 2012J), with photo-x; tests presented in 


Hildebrandt et al 


iBeniamin et al J d2013l) . Details on the CFHTLenS analysis pipeline 
and cosmology-independent systematic tests are presented in 
iHevmans et all 120121) . 

To subdivide the foreground lens galaxies in our analy- 
sis we employ stella r mass estimates c omputed using LePhare 
dArnouts et all 1 19991 : lllbert et all l2006ll . They are based on the 
stellar population synthesis (SPS) package of iBruzual & Charlot 
1 2003) and assume a stellar initial mass function from Ichabriei 
1 20031) . From this we built 18 templates with two different metal- 
licities and nine exponentially decreasing star formation rates. We 
allow the du st ex tinc tion to vary between 0.05 to 0.3 mag us¬ 
ing a ICalzetti et al J d2000l ) extinction law, and 57 starburst ages 
ranging from 0.01 t o 13.5 Gyr (for more details see Sec. 2.1 of 
IVelander et al.ll2014h . 


3.2 KSB shapes for bright foreground galaxies 


The lensfit shape measurement algorithm has been optimised to 
obtain accurate shape estimates for the typicall y fain t and_only 
moderately resolved distant source galaxies dMiller et alj 120131) . 
as n eeded for an unbiased cosmological weak l ensing analysis 
(e. g. iKi lbinger et~ah 1 120131: IHevmans et al] 120131; iFu et al.l I20l4 
Kitching et al . 12014 ). For our analysis of galaxy halo shapes we ad¬ 
ditionally require ellipticity estimates for the foreground lens galax¬ 
ies. Many of these are fairly bright and extended, which may result 
in an exclusion of a galaxy during the lensfit shape analysis. This 
can be caused either because of a size comparable or exceeding the 


6 For example, in our analysis of haloes from the Millennium Simu¬ 
lation (see Sect. [4} that are populated by early type galaxies, we find 
(l e h|) = 0.163, and only 0.5% of the haloes have zy, > 0.5. 
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employed postage stamp size (48 pixels), or outer isophotes signif¬ 
icantly overlapping with neighbouring galaxies, or the presence of 
substantial substructure in the galaxy, which is not well described 
by the employed bulge+disk model (see lMiller et al.|[2013l) . 

In order to not exclude these galaxies from our sample of lens 
galaxies, we perform shape measu rements using th e KSB algo¬ 
rith m originally propos ed bv lKaiser. Squires & Broadhursll ( 1995 ) 
and lLuppino & Kaiser! J 19971 ). Based on weighted brightness mo¬ 
ments this shape measurement algorithm is more robust to the pres¬ 
ence of resolved substructure or nearby galaxies. Here we e m¬ 
ploy th e i mplementation of the algorithm deta iled in lHoekstra et alj 
(119981) and lHoekstra, Franx & Kuiikenl d2000h. which was t ested in 
the blind challenge s of the STEP project IHevmans et alj|2006al ; 
iMassev et al.l 120071) and also employed for earlie r weak lens- 
ing analyses of CFHTLS data dHoekstra et all 1 20061) . While there 
are indications for remaining residual systematics in earlier cos- 
mological weak lens i ng studies of CFHT LS data using KSB 
dKilbinger et alj 120091 : IHevmans et aid 120121) . we note that these 
are expected to mostly originate from poorly resolved, low signal- 
to-noise galaxies. In contrast, the studied lens galaxies have high 
signal-to-nois^H and are w ell-resolved. Thus, they are less sensitive 


to noise-related biases (e.g. 

Melchior & Viola 2012; Refregier et al. 

2012 ; ; Kacprzak et al. 2012 

; Viola, Kitching & Joachimi 2014) and 

require smaller PSF corrections. Also see Sifon et al. ((2015!) 


where the KSB results f or bright galaxies are compared to GAL- 
FIT shapes dPeng et ahl 12002 ). Furthermore, as demonstrated in 
Sect. 12.21 additional ellipticity correlations between lenses and 
sources due to imperfect PSF anisotropy correction or cosmic shear 
cancel out in the analysis at the relevant scales thanks to the em¬ 
ployed estimator. Thus, the application of KSB shapes for lens 
galaxies without lensfit shapes does not compromise the system¬ 
atic accuracy of our measurement, but only increases the statistical 
precision. Note that differing from earlier KSB studies of CFHTLS 
data we conduct galaxy shape measurements on individual expo¬ 
sures and not stacks, similarly to the lensfit analysis, and combine 
the shape estimates on the catalogue level as the weighted mean 
estimate. 


3.3 Lens sample 

To obtain a sample of foreground lens galaxies we pre¬ 
select relatively bright objects (i < 23.5) that are resolved 
(CLASS.STAR < 0.5, star_flag = 0, non-zero shape weights 
from lensfit or KSB) and feature a single-peaked photomet¬ 
ric redshift probability dis tribution function (ODDS > 0.9, 
see IHildebrandt et al.112012 ). We select lenses within the range 
of best-fitting photometric redshifts 0.2 < Zh < 0.6, which we 
subdivide into four lens redshift slices of width Azb = 0.1. 
We split the galaxies into red (Tbpz ^ 1.5) and blue lenses 
(1.5 < Tbpz < 3.95) using the photometric type Tbpz from BPZ. 
In order to approximately sort the lenses by halo mass and op¬ 
timise the measurement signal-to-noise ratio, we also subdivide 
them according to stellar mass log 10 M, as detailed in Table |T| 
When measuring the anisotropic shear signal, contributions from 
different lenses are weighted according to the lens ellipticity (see 


' In the highest-redshift redshift slice (0.5 < z < 0.6), the median signal- 
to-noise ratio defined as FLUX_AUTO/FLUXERR_AUTO from SExtractor 
is 181 (193) for the blue (red) lenses in the lowest stellar mass bins consid¬ 
ered. 


e.g. Eq. [J5J. Here we restrict the analysis to lenses in the well- 
constrained ellipticity range 0.05 < |e g | < 0.95. 

An important aspect of our study is the direct comparison 
of our measurements from CFHTLenS to results from mock data 
based on the Millennium Simulation (see Sect.0- For the mock 
data we only have measurements for central haloes but not for 
satellites (see Sect. 1421 To ensure that the results are compara¬ 
ble, we aim to minimise the fraction of satellites in our CFHTLenS 
lens samples. In addition, our model assumes that the anisotropic 
shear signal is caused by isolated elliptical NFW mass distributions 
(Sect.[2}. This may be a reasonable approximation for centrals, but 
it is likely a poor descriptions for satellites embedded into a larger 
halo._ _ 

IVelander et al.l d2014l) fit the isotropic galaxy-galaxy weak 
lensing signal around red and blue lenses in CFHTLenS using a 
halo model approach. They find that the satellite fraction is typ¬ 
ically low both for blue lenses and for red lenses which have a 
high stellar mass. However, it increases steeply for red lenses to¬ 
wards lower stellar mass. To reduce the fraction of satellites in our 
sample, we therefore generally exclude red lenses with stellar mass 
log 10 M* < 10. For the lowest stellar mass bin included in our 
analysis for red lenses (10 < log 10 AT* < 10.5), IVelander et alj 
( I 2 OI 4 I) estimate a satellite fraction of a = 0.23 ± 0.02. To fur¬ 
ther reduce the fraction of satellites, we apply an additional cut 
to the galax ies in this bin usi ng the internal flag from SEXTRAC¬ 
TOR ( iBertin & Arnoutslll996l) . This removes any galaxy which is 
either partially blended with another object, or which has a nearby 
neighbour possibly affecting the measurement of the MAG_AUTO 
magnitude. Many of these galaxies are located close to a brighter 
early type galaxy, as expected for satellites. We do not filter on this 
flag for the other stellar mass bins. In particular, this would exclude 
many bright early type galaxies, which are presumably centrals but 
have faint nearby neighbours. 

Given that the remaining fraction of satellites in our lens sam¬ 
ple is low, we expect that they have a negligible impact on our 
result^. As consistency check for this we investigate the impact 
of the lens bin with the highest satellite fraction on our joint con¬ 
straints in Sect. l3.5.2l 


3.4 Source sample 

From all galaxies with lensfit shape measurements and non¬ 
vanishing shape weights, we select those with a best-fitting pho¬ 
tometric redshift in the range zi iUp per + 0.1 < Zb < 1.3 as our 
source background sample, where zi, upper indicates the upper limit 
of the corresponding lens redshift slice. Above the upper limit at 
Zb = 1. 3 the uqriz CFHTLenS photometric redshifts become un¬ 
reliable ( IHildebrandt et aljr2012l ; lBeniamin et alj|2013l) with large 
redshift uncertainties and a partial contamination by low-redshift 
galaxies. Thus, we exclude these galaxies from our analysis to be 


8 A net alignment of satellites with respect to the tangential shear field 
of their parent halo could introduce spurious signal. However, we expect 
that this has a negligible impact on our study for the followi ng reasons: 
First, t he satellite fraction is low for our lens samples. Second, Sifon et al. 
<2015h place tight limits on the net alignment of cluster galaxies, showing 
that it must be very weak. Given our selection of lenses with relatively high 
stellar mass, we expect that many of the satellites present in the sample are 
actually cluster members. Finally, our analysis using (/ — / 4 s)AE reduces 
the impact of spurious alignment between lenses and sources, which also 
applies here to first order. 
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Table 1. The lens galaxy samples. 


3.5 Results 


Colour 

Stellar mass [Mq] 

N( 0.2 s; z, < 0.6) 

Ce 

Red 

10 < log 10 M* < 10.5 

81763 

0.35 

Red 

10.5 < log 10 M* < 11 

93032 

0.30 

Red 

log 10 M * > 11 

25982 

0.23 

Blue 

9.5 < log 10 M* < 10 

166604 

0.38 

Blue 

10 < log 10 M » < 10.5 

91692 

0.36 

Blue 

log 10 M* > 10.5 

33612 

0.30 


Note. — Overview over the sub-sample of lens galaxies used: Column 1: 
Split between red (Tbpz C 1-5) and blue lenses (1.5 < Tppz < 3.95) 
using the photometric type Tfipz front BPZ. Column 2: Stellar mass 
range. Column 3: Number of selected lenses in the redshift interval 
0.2 ^ z\ < 0.6. Column 4: Ellipticity dispersion of the selected lenses 
with 0.05 < |e g | < 0.95 combining both ellipticity components. 


conservative. As done by IVelander et al j l l2014h . we further opti¬ 
mise the separation of lenses and sources by also removing source 
galaxies whose 95% redshift confidence regions computed by BPZ 
overlap with the lens redshift slice. When we compute the galaxy- 
galaxy lensing signal we weight the contributions from the individ¬ 
ual source galaxies according to their effective geometric lensing 
efficiency /3f s = J /3(zi, c , z s )pi(z a )dz s . This is estimated from the 
full photo metric redshift probabil ity distributions of the sources 
Pi(Zs) (see lHildebrandt et al.l2012h and uses the central of the thin 
lens redshift slice z\, c . 

Miller et al.| j2013l) study the impact of noise bias for lensfti 
shape measurements and derive an empirical correction for multi¬ 
plicative bia s as a function of ga laxy size and signal-to-noise ratio. 
Following IVelander et alj ( l2014h we account for the net effect of 
this bias on the estimated shear profiles taking the individual shape 
weights and lensing efficiencies of the sources into account. We 
note that this bias correction cancels out for the halo ellipticity con¬ 
straints except for its impact on the estimate of r 2 ooc, which is used 
as upper limit for the fit range and to estimate r s for the anisotropic 
shear field model. _ 

iHevmans et alj d2012l) conduct a number of non-cosmological 
tests to flag CFFITLS fields which are likely affected by residual 
shape systematics. To be conservative, we base our primary analy¬ 
sis on the 129 out of 171 fields which pass these systematics testa 9 10 ! 
(“pass fields”), but for comparison we also provide constraints ob¬ 
tained from all fields. As explained in Sect. 12.11 the impact of both 
residual shape measurement systematics and foreground cosmic 
shear cancel out in our analysis as long as their shear correlation 
function £_ is sufficiently small on the scales of interest. As we will 
see in Sect. l4.3.4l this approach leaves negligible residuals from the 
expected cosmic shear signal. Accordingly, we expect that the fields 
classified as being useful for cosmic shear measurements can also 
be used to derive robust constraints on the aligned halo ellipticity. 


9 We use the centre of the thin lens redshift slice instead of the individual 
best-fit lens redshift for computational efficiency. We note that the influence 
is negligible for our study as verified by varying the width of the slices. 

10 See http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/ 

community/CFHTLens/README_catalogs_release.txt. 


3.5.1 Measured signal and fit range 


For each set of lens galaxies we measure the isotropic and 
anisotropic shear signal in 25 logarithmic bins of transverse phys¬ 
ical separation between 20kpc//i7o and 1.2Mpc//i7o. We compute 
the signal separately for the lenses in each ~ ldeg 2 CFFITLS field, 
where we however use larger cut-outs from the mosaic catalogue 
of source galaxies to ensure a good coverage for the lenses close 
to the edge of a CFFITLS field. We then combine the signal from 
all fields according to their weight sums in each radial bin. Simi¬ 
larly, we compute a combined signal from all lens redshift slices 
for each field and the full survey. Figs.[3] and [4] show the signal 
combining all redshift slices for the different stellar mass bins for 
red and blue lenses, respectively. Error-bars are computed by boot¬ 
strapping the CFFITLS fields contributing to the combined sur¬ 
vey signal. In this study we are primarily interested in the rela¬ 
tive asymmetry in the shear field. For detailed investigations of 
the isotropic galaxy-galaxy lensing signal in CFH TLenS please 
Velander et al. 11201 4l) : [Hudson et alj ( 120151) . Coupon et alj 


d2015l) . 

Before we can obtain constraints on the halo ellipticity we 
have to choose the radial range which is included in the fit. We are 
interested in the shape of the dark matter-dominated matter halo 
surrounding the lens. Thus, we exclude small scales which are ex¬ 
pected to be affected by the baryonic component of the lens. This 
is visible from the isotropic galaxy-galaxy lensing signal (see the 
top panels of Figs.[3]and[4]l, showing a strong excess signal at small 
scal es compared to the NFW fit constrained fr om larger scales (see 
also lVelander et al1l2014l : lHudson et alj|2015l who include a bary¬ 
onic component in the fit to the isotropic signal). 

In addition, the signal at very small scales might be influenced 
by systematic effects in the source detection and shape measure¬ 
ment process originating from the presence of the nearby bright 
foreground lens galaxies. As a first check for this we studied the rel¬ 
ative azimuthal variation in the source density. Fig.[5]shows the ra¬ 
dial dependence of YLiLi cos (2A0;)//V for all source-lens pairs in 
the highest stellar mass bins for both the blue and red lens samples, 
combining all redshift slices. The decrement at r < 45kpc//i7o 
shows that we have a higher source density in the direction of 
the lens minor axis than the lens major axis. This suggests that 
either the object detection and deblending, or the shape measure¬ 
ment is more effective in the direction where the lens light has a 
smaller impact, which is not surprising. To ensure that this cannot 
affect our measurements we therefore only include radial scales 
r > 45kpc//i7o for all lens bins, which also matches our goal to 
remove the scales where baryons are important. In Appendix lAlwe 
present an additional test for shape measurements close to bright 
lens galaxies using simulated galaxy images. It shows that any re¬ 
sulting spurious signal should be small compared to the statistical 
uncertainties from CFHTLenS within our fit range. 

We also have to select an outer radius for the radial range 
which we include in the halo shape analysis. Here we decide to 
include the signal out to V 2 ooc. At significantly larger radii the 
isotropic galaxy-galaxy lensing signal shows an excess signal for 
the lower stellar mass bins (see the middle and right panels in the 
top row of Fig.[3] showing exce ss sig nal at r > 1.6r 2ooc). Inter¬ 
preted in the halo model (see e.e. IVelander et al] 2014h . this excess 
signal is a combination of the satellite term caused by the central 
galaxy (for lenses which are satellites) and the two-halo term from 
neighbouring haloes. Hence, at these scales the halo of the lens¬ 
ing galaxy can no longer be regarded as isolated, which would re- 
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Figure 3. Measured isotropic (top row of panels) and anisotropic (rows two to four, note the varying ?/-axis scale) shear signal around red lenses in the 
CFHTLenS fields passing the systematics tests for cosmic shear as function of radial distance r. The anisotropic signal has been scaled by r for better readabil¬ 
ity, where rows two, three, and four show the signal components (/ — /isjAE, /AE, and — AS, respectively. From left to right, we show the stellar mass 
bins log 10 M« > 11, 10.5 < log 10 M* < 11, and 10 < log 1() M* < 10.5, combining all lens redshift slices. For the isotropic signal the curve shows the 
best-fitting NFW shear profile constrained within 45 kpc//? 7 o < r < 200 kpc//i 7 o. For (/ — / 4 s) AS the curves show models corresponding to the best-fit 
isotropic model and the best-fitting j), (solid curves), as well as /), £ {+1,0, —1} (dotted curves) for comparison. For /AS and —/ 4 s AS, the dashed 
curves show model predictions for the best-fitting /n- The best-fitting /), has been constrained from (/ — // 15 ) AS and AS within 45 kpc /hjo < r < r 200 c 
(indicated by the vertical dotted lines), with r 2 oo c estimated from the fit to the isotropic signal. 


quire a much more complicated modelling scheme. We note that 
the galaxy bins with higher stellar mass show only a weak excess 
signal at large radii, presumably due to a low satellite fraction. So 
while we could in principle extend the fit range for these galaxies 
slightly, we decided to use r 2 ooc as upper limit for all bins to keep 
the analysis more homogeneous. 

Our approach requires that we first obtain an estimate for 
f 200 c- For this we fit the isotropic signal in the fixed radial range 
45 kpc//i 7 o < r < 200 kpc/fi. 70 . Here we use the central redshift 
of our lenses z\ = 0.4 when computing the isotropic NFW model 
for the combined analysis of all redshift sliceiTH. 

For comparison, we also plot the components of the 
anisotropic signal /AE and —/ 45 AE in the bottom panels of 
Figs.[3] and [4] We note that the blue galaxies shown in Fig.[4]show 
a tendency for /AE < 0 and —/ 45 AE > 0, especially towards 


11 The central lens redshift agrees well with the effective lens redshift. The 
latter is computed according to the weight sums of the thin redshift slices 
for a fixed radial bin, yielding effective lens redshifts of 0.38 to 0.41 for 
the different galaxy types and stellar mass bins. We verified that the exact 
choice of the model z\ has minimal impact on our j\ x constraints. 


larger radii, which is consistent with the expected trend for cosmic 
shear. We will discuss this further in Sect.[5] 

From the bootstrapping analysis we find that off-diagonal 
terms in the correlation matrix are small, justifying our analysis 
approach (see Sect. l2.1.3l . Within the fit range of the isotropic sig¬ 
nal the average of the off-diagonal elements is consistent with zero 
at the ~ 1 — 2 cr level for all lens bins with (cor; i ,);> ; ,j < 2 %. 

3.5.2 Constraints on fh 

The results of the fits to the CFHTLenS data are presented in Table 
[2]and we show the model fits to the data in Figs.[3]and[4] For each 
lens colour and stellar mass bin we fit the combined shear signal 
from all redshift slices to ensure that the isotropic shear profile, 
which determines r 2 ooc, is measured with high significanca 12 !. 

For none of the individual lens bins do we detect an /h sig¬ 
nificantly different from zero. We also compute joint constraints 

12 For comparison we repeated the measurement where we initially anal¬ 
yse each redshift slice separately and combine the constraints when esti¬ 
mating /'|,. This led to nearly identical, and within the statistical uncertainty 
fully consistent results. 
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Figure 4. Measured isotropic (top panels) and anisotropic (rows two to four, note the varying y-axis scale) shear signal around blue lenses in the CFHTLenS 
fields passing the systematics tests for cosmic shear as function of radial distance r. From left to right, we show the stellar mass bins log 10 M* > 10.5, 
10 < log 10 M* < 10.5, and 9.5 < log 10 -&/* < 10, combining all lens redshift slices. For further details see the caption of Fig.[3 


Table 2. Weak lensing results using the CFHTLenS data. 


Colour 

Stellar mass 

r* 200 c 

37200c 

pass 

fields 

all fields 


[M e ] 

[kpc//i 70 ] 

[10 11 M Q /h 70 \ 

/h 

X 2 /d.o.f. 

/h 

X 2 /d.o.f. 

Red 

10 < log 10 M, < 10.5 

170 

8.6 ± 1.7 

—0 11+°- 73 
U-il —0.73 

9.6/7 

—0 3i+ 0 - 59 
-0.60 

7.8/7 

Red 

10.5 < log 10 M, < 11 

253 

28.2 ±2.9 

0 f)i+°- 35 
U-U± -0.35 

9.5/10 

—0 20" 1 ” 0 ' 29 
U.ZU-O 29 

8.8/10 

Red 

log 10 M * > 11 

377 

93.3 ± 10.3 

—0 nq+°- 38 
u - uy -0.38 

11.3/12 

—0 ni+°- 33 
U ‘ U1 -0.33 

14.0/12 

Red 

log 10 M* > 10 



U.U^_o.25 

30.5/31 

-0.17±g;li 

30.6/31 

Blue 

9.5 < log 10 M* < 10 

128 

3.7 ±0.8 

-i 

-0.81 

7.5/5 

0.89l°j® 

11.2/6 

Blue 

10 < log 10 M, < 10.5 

191 

12.1 ±2.0 

U - iy -0.50 

2.8/8 

0 35+ 0 - 49 
U.OO-0.49 

3.8/8 

Blue 

log 10 M* > 10.5 

234 

22.2 ±4.4 

1 11+0.67 

r.11—o.64 

7.3/9 

0 51 +°- 52 

U.C)l_o.5i 

9.9/9 

Blue 

log 10 M» > 9.5 



0 69+ 0,37 

u.oy _ 0 36 

20.5/24 

0 '56+ 0 - 34 
u - OD -0.33 

24.8/25 


Note. — Summary of the results from the analysis of the CFHTLenS data: Column 1: Split between red and blue galaxies. Column 2: Stellar mass range. 
Columns 3 and 4: Radius r 200 c and mass M 20 OC as estimated from the isotropic component of the shear signal from the fields passing the systematics tests 
(“pass fields”). Columns 5 and 6: Aligned ellipticity ratio and reduced x 2 for the fit from the analysis of the “pass fields”. Columns 7 and 8: Aligned 
ellipticity ratio /h and reduced x 2 for the /h fit from the analysis of all fields. 


from the (/ — f± 5 ) AE and AE profiles of all stellar mass bins as 
explained in Sect. l2.1.3l yielding /h = —0.04 zb 0.25 for the red 
lenses and /h = 0.69^o 36 for the blue lenses when restricting the 
analysis to the 1 29 fie lds passing the systematics tests described in 
iHevmans et alJ <120121 . “pass fields”). For comparison, we estimate 
/h = —0.17 ± 0.21 for the red lenses and /h = 0.561q 33 for the 
blue lenses when including all 171 CFHTLenS fields. In Table [2] 
we also list reduced \ 2 values which suggest that the models fit the 


data reasonably well in both cases when considering all lens bins 
together, but we note slightly lower x' 2 /d.o.f. for the blue lenses 
when using the “pass fields” only. 

As a consistency check for the possible impact of satel¬ 
lite galaxies in the lens sample (see Sect. l3.3b we also com¬ 
pute joint constraints for the red lenses now excluding the 
stellar mass bin with the highest expected satellite fraction 
(10 < log 10 M* < 10.5). In this case we obtain almost un- 
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Figure 5. Measurement of the relative azimuthal variation in the source den¬ 
sity in CFHTLenS with respect to the lens major axis as function of radial 
distance. The open red circles (blue crosses) show JT =1 cos (2A9i)/N 
for the red (blue) lenses in the correspondingly highest stellar mass bin com¬ 
bining all lens redshift slices. The decrement at r < 45 kpc/hyo indicates 
a higher source density in the direction of the lens minor axis compared to 
the major axis. To ensure that this cannot influence our analysis we only 
include scales r > 45 kpc/hzo in the halo shape analysis. 


changed results /h = —0.02 ± 0.26 for the “pass fields” and 
/h = —0.14 ±0.22 for all fields, suggesting that the remaining 
satellites have a negligible impact on the joint constraints. 


4 ANALYSIS OF SIMULATED DATA BASED ON T HE 
MILLENNIUM SIMULATION 

To better understand the halo shape signal that we should expect 
from non-idealised haloes and in the presence of misalignments 
between galaxies and their matter haloes, we analyse a simulated 
data set based on ray -tr acing through the M illennium Simulation 
llSpringel et a 1.1 [200511 bv IHilbert et alJ d2009h . This also allows us 
to test the correction for systematic shear (see Sect. l2.lt in the 
presence of a realistic cosmic shear field. We refer the reader to 
iHilbert et alj (2009 ) regarding the details of the ray-tracing. Here 
we only summarise some of the main characteristics relevant for 
our analysis. 


4.1 Mock shear catalogues 

The simulated catalogues comprise 64 light cones, each with 
an area of 4 x 4deg 2 , which we subdivide into patches of 
1 deg 2 to facilitate a bootstrap analysis similar to the CFHTLenS 
fields. We use all galaxies at 0.65 < z < 2.15 as source sam¬ 
ple (z m edian = 1.187), providing a high source density of 
49.6/arcmin 2 . We wish to obtain high S/N estimates of the sim¬ 
ulated halo shape signal and therefore add only small shape noise 
with g e = 0.03, which was chosen to be of the same order as the 
noise introduced by cosmic shear. This low level of shape noise 
leads to stronger relative noise contributions from cosmic shear 


compared to the CFHTLenS analysis. This increases the noise cor¬ 
relations between different radial bins. While they can be sub¬ 
stantial for the isotropic signal AE for some of the lens bins 
with 5% < |(corjj)j>_,j < 45%, they are generally small for the 
anisotropic signal (/ — / 4 s)AE with |(cori 1 j)i>j| < 4%. Noise 
in the latter dominates the uncertainties in the /h constraints. We 
therefore expect that the net impact of the noise correlations on 
our constraints is negligible. As a consistency check we also repeat 
the analysis of the simulated data with more realistic shape noise 
(<r e = 0.25). Here we find consistent results, but the increased sta¬ 
tistical uncertainty makes it impossible to detect the signal for some 
of the weakly aligned lens models. 


4.2 Mock lens galaxies 


For the foreground len s es we make us e of g alax y shapes com¬ 
puted b y ljoachimi et al] d2013al her eafter J13a) and Joachimi et alj 
d2013bL hereafter |j!3bl) . Following iBett et a" d2007h . haloes and 
their member candidate part icles are fi rst identified by a friends- 
of-friends algorithm foavis et al.ll 19851) . In a second step, merger- 
tree data is used to remove particles belonging to substructures that 
are only temporarily in the halo vicinity. Halo shapes are then es¬ 
timated via the full quadrupole tensor of the halo’s matter distri¬ 
bution that remains after removal of these transients. This infor- 
mat ion is comp lemented with semi-analytic galaxy evolution mod¬ 
els l iBower et al.1120061) and a classification of galaxy morphologies 
into early and late ty pes via the bulge-to- total r atio of the rest-frame 
A'-band luminosity JParrv. Eke & Frenkll2009t) . They also separate 
galaxies into centrals (galaxies in the most massive substructure of 
a halo) and satellites. Here we only use f oreground galaxies clas¬ 
sified as centrals, because |j!3al and |j!3bl had to assign simplistic 
galaxy shapes to the satellites given that their dark matter haloes 
were poorly resolved. Also, we only keep those foreground galax¬ 
ies with a sufficient particle number for robust shape estimation, 
see lj 1 3al For the centrals J13alan d ij 13bi assign shapes adopting the 
scheme of lHevmans et all i2006bl) . Here, late type disc-dominated 
galaxies are aligned such that their spin vector is parallel to the 
angular momentum vector of their host dark matter halo. These 
shapes are then projected onto the plane of the sky. For early type 
galaxies it is assumed that the shapes of the galaxies follow the 
shapes of their dark matter haloes. For this. lJ 13al project the ellip¬ 
soid given by the eigenvectors and eigenvalues of the inertia tensor 
for each halo on to the plane of the sky and use the resulting ellipse 
as the shape of the galaxy. 

We additionally consider galaxy ellipticities that are mis¬ 
aligned with respect to their host dark matter halo as detailed in 
IJ13a . For early type galaxies the misalignment angles were drawn 
fro m a Gaussian distribution with a scatter of 35° as estimated 
by lOkumura. Jine &~Lil d2009l ) from the distribution of satellites 
around luminous red galaxies (LRGs) in SDSS. For late type galax¬ 
ies Ijl3bl emplo y a m isalignment distribution based on a fitting func¬ 
tion that lBetd 120121) determined using a compilation of_simulated 
haloes with baryons and galaxy formation physics I Bett et a l.l20ld : 
iDeason et all201 it ICrain et al]|2009l : lOkamoto et al,ll2003) . 

Similarly to the analysis of the CFHTLenS data, we select lens 
galaxies in the redshift range 0.2 < z\ < 0.6 split into two redshift 
slices (which we analyse separately because of better S/N com¬ 
pared to CFHTLenS) and stellar mass bins, see Table[3]for details. 
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12 Schrabback et al. 


Table 3. Weak lensing results using the mock data based on the Millennium Simulation, both for the analysis with and without foreground cosmic shear (c.s.) 
applied. 


Type 

Stellar mass 


^1 , max 

o- e 

r 200 c 

37200c 

fh (aligned) 

f h (misaligned) 






[kpc//i 73 ] 

[lO U Jlf 0 /ll73] 

with c.s. 

without c.s. 

with c.s. 

without c.s. 

Early 

9.5 < log 10 M, < 10 

0.2 

0.4 

0.19 

168 

7.84 ± 0.03 

0 504 4 ” 9 027 

u.ou ^_ 0 025 

0 511 4-0,026 
U.O1-L—0.026 

0.2031°;°!! 

0.2191°;°!° 

Early 

10 < log 10 M, < 10.5 

0.2 

0.4 

0.19 

226 

18.97 ± 0.06 

0 5QC+0.015 

U.DoD_o Q i5 

0 644"*" 0 015 
u.on±t_ 0 Q15 

0.3061°;°!° 

0-2981°;°!° 

Early 

10.5 < log 10 M» < 11 

0.2 

0.4 

0.19 

375 

86.19 ±0.25 

0.759t°;“! 

fl 779+ 0 - 010 
u -' ‘ z -0.008 

0 Q 7 Q + 0-010 

u.o i o _ 0 008 

0 370 4-0 010 
u.oi u_o 008 

Early 

11 < log 10 M» < 11.5 

0.2 

0.4 

0.21 

563 

291.14 ± 1.69 

0.846±°;°i? 

0.8631°;°!! 

0 402”*” 0 011 
U.^UZ_o 010 

0 407”*”°'° 12 
u -^ u ‘- 0.010 

Early 

9.5 < log 10 M* < 11.5 

0.2 

0.4 

0.19 



0 71 5 4-0 - 007 
u. i ^o_ 0> 006 

0.726l°;°Qg 

0 Q44 + 0 007 
u - c> ^-0.006 

0 340 4-0 007 
U.O^O _ 0 Q06 

Early 

9.5 < log 10 M, < 10 

0.4 

0.6 

0.20 

156 

7.56 ± 0.03 

0 407”*” 0 ' 021 
u -^ u ' - 0.021 

0 463+ 0 022 

u.huo _ 0 020 

0 144 4 " 0 022 

U-J-^- 0.021 

0 209”*” 0 ' 022 
u.zuy _ 0 02 i 

Early 

10 < log 10 M, < 10.5 

0.4 

0.6 

0.20 

213 

19.15 ± 0.05 

0 501+ 0 - 013 
U.OUl_o 012 

0 535 + 0 013 

u.ooo_o 012 

n 909 + 0.013 
u.^o^_ 0.012 

0 259”*” 0 012 

u.zoy _ 0 012 

Early 

10.5 < log 10 M» < 11 

0.4 

0.6 

0.21 

351 

86.34 ±0.24 

0-660l°;“l 

0.6791°;°°° 

0 29Q 4 " 0 009 

u.zyy _ 0 007 

n 009 + 0.009 

U.ooz _ 0 007 

Early 

11 < log 10 M» < 11.5 

0.4 

0.6 

0.22 

502 

251.61 ± 1.37 

0-8051°°!;; 

0.8271°;°!! 

0 35Q + 0 011 

u.oo »_ 0 oio 

0 37Q 4-0 011 
u.o/y _ 0 010 

Early 

9.5 < log 10 M* < 11.5 

0.4 

0.6 

0.21 



0 oo7+° °06 

U.OU i —0.006 

0 «Q «+ 0 006 

U.O 0 O_o 006 

fl 979 + 0 - 007 

u.ziz_Q 005 

0 qqo 4-0 006 
u.ouo _ 0 005 

Early 

9.5 < log 10 M, < 11 

0.2 

0.6 

0.20 



0 61 fi 4-0 - 006 

U.OlO_0 005 

0.6371°;°°° 

0 OS'S ”*" 0 006 

u.zoo _ 0 004 

0 QQ7 4-0 - 005 
U.OU/_o 005 

Early 

9.5 < log 10 M* < 11.5 

0.2 

0.6 

0.20 



0 o57 "*" 0,005 

u.uo/ _ 0 Q04 

0.6781°;°°! 

0 QQ4 + 0 005 

U.oUlt_o 004 

0 324" 1 ” 0 ’ 004 

u.oz ^_ 0 004 

Late 

9.5 < log 10 M» < 10 

0.2 

0.4 

0.32 

131 

3.68 ± 0.01 

0 073”*” 0,016 

U.lUO_o 016 

0 . 0681 °;°!! 

0.0341°;°!! 

0 026”*” 0 ’ 017 
u.uzo _ 0 016 

Late 

10 < log 10 M» < 10.5 

0.2 

0.4 

0.32 

171 

8.17 ± 0.03 

0 -Jos -1 " 0 - 014 

u.izo_o Q13 

n i qq+o.ois 
u. -Loo— 0,013 

0.0381°;°!! 

0.0641°;°!! 

Late 

10.5 < log 10 M, < 11 

0.2 

0.4 

0.32 

245 

24.02 ± 0.10 

0.1291°;°!! 

0 1 qo+ 0.012 

u.±oz _ 0 on 

0 0 S 1 4-0 012 
U.UOl _0 012 

0.0591°°!! 

Late 

9.5 < log 10 M* < 11 

0.2 

0.4 

0.32 



0 . 1081 °;°°° 

O.llllaoS! 

0 QQ 4 4-0 009 
U.UO4 _ 0 008 

0 049”*" 0 009 
u.u^y _ 0 008 

Late 

9.5 < log 10 M, < 10 

0.4 

0.6 

0.32 

124 

3.82 ± 0.01 

0 041 4-0,014 
u - u ^ :i - 0.012 

0 . 0681 °;°!° 

0 007”* -0 ' 014 
U.UUf _ 0.012 

0 013 4-0 - 014 
u.uio _ 0 013 

Late 

10 < log 10 M, < 10.5 

0.4 

0.6 

0.32 

163 

8.68 ± 0.03 

0 101 4-0012 
u.iui_o oio 

n i OO+ 0.012 

U.±zz _ 0 on 

0 022”*” 0 ' 012 
u.uzz _ 0 on 

0 QQ 1 4-0 012 
u.uoi_o.on 

Late 

10.5 < log 10 M» < 11 

0.4 

0.6 

0.32 

230 

24.38 ± 0.08 

0-1151°;°!° 

n 19O+0-010 
U.lzo_o oio 

0 029”*" 0-011 
u.uzy _ 0 oio 

0 Q 37 4-0 011 
u.uo / —0.009 

Late 

9.5 < log 10 M* < 11 

0.4 

0.6 

0.32 



0.085l°;°°! 

0 1Q4 + 0 008 
u - lu ^-0.006 

0 01 q 4-0 007 
u.uiy_ooo 6 

0 Q27 4-0 - 007 
U.UZ 1 — 0.007 

Late 

9.5 < log 10 M* < 11 

0.2 

0.6 

0.32 



0.095l°;°°! 

o 

o 

1 + 

o o 

o o 

o o 

05 

0 Q25 4-0 006 

u.uzo _ 0 004 

0 OS6 4-0 ' 006 
U.UOD_ 0> 005 


Note. — Summary of the results from the analysis of the Millennium Simulation. Column 1: Split between early and late types. Column 2: Stellar mass 
range. Column 3: Minimum lens redshift. Column 4: Maximum lens redshift. Column 5: Ellipticity dispersion of the selected lenses with 
0.05 < |e g | < 0.95 combining both components. Columns 6 and 7: Radius r 200 c and mass M 20 OC as estimated from the isotropic component of the shear 
signal for the aligned lens models with cosmic shear (very small differences for the other models). Columns 8 and 9: Aligned ellipticity ratio /[, for aligned 
lens models (models Est and Sal from |jl35 ) with and without foreground cosmic shear applied, respectively. Columns 10 and 11: Aligned ellipticity ratio 
J\ ] for misaligned lens models (models Ema and Sma Irom l.l 13 fal l with and without foreground cosmic shear applied, respectively. 


4.3 Results 

We summarise the main results of the analysis based on the Mil¬ 
lennium Simulation in Table [3] listing estimated halo masses, as 
well as /h constraints for both aligned and misaligned galaxy shape 
models, each with and without foreground cosmic shear applied. 
For this analysis we use modified ACDM cosmological parameters 
for the evaluation of angular diameter distances, matching the in¬ 
put to the Millennium Simulation, with = 0.25, 12 a = 0.75, 
Ho = 73/173 km/s/Mpc. 


4.3.1 Halo masses 


In the CFHTLenS analysis we use stellar mass as proxy for halo 
mass. For the analysis of the Millennium Simulation we again 
split the lenses into stellar mass bins. However, these stellar mass 
estimates are base d on semi-analytic galaxy evoluti on models 
{Bower et al. 200(J) which have some uncertainty (e.g. iKim et all 
120091: ISanhiba et allhoi 21) . For a direct comparison of the results 
from CFHTLenS and the simulation it is more relevant that the 
included bins roughly span the same range in halo mass. Here 
we note that the halo mass estimates inferred from the isotropic 
shear signal agree reasonably well for the late type (blue) lenses. 


whereas the early type (red) lenses yield higher halo masses in 
a given stellar mass bin in the simulation. To roughly match the 
halo mass range when computing joint constraints on fh from 
the simulation, we therefore remove the highest stellar mass bin 
(11 < log 10 M* < 11.5) and instead include one additional, lower 
stellar mass bin (9.5 < log 10 M* < 10). Given the only moderate 
dependence of /h on mass (Table[3]i- this approximate matching is 
fully sufficient for our goal to provide model predictions. 


4.3.2 Shear profile plots 

We plot the isotropic and anisotropic shear profiles for an illus¬ 
trative subset of the combinations of lens bins and shape mod¬ 
els in Figs.[6] and |7] From top to bottom the panels show the 
isotropic profile AE(r) and the corresponding anisotropic profiles 
(/ — /isjAEjr), /AE(r), and —/ 45 AE(r), both with (crosses) 
and without (circles) foreground cosmic shear applied to the lens 
ellipticitieEl 

Fig.[6] shows early type lenses with high stellar mass 


13 The foreground cosmic shear is always applied to the sources, but it only 
has a net effect if it is applied to both lenses and sources. 


© 2015 RAS, MNRAS 000. [111211 




















CFHTLenS: Halo ellipticity constraints 13 


10.5 < log 10 M* < 11. Here, the left and middle columns show 
the signal of lenses with aligned shape models at lower 
(0.2 < z\ < 0.4) and higher (0.4 < z\ < 0.6) redshift, respec¬ 
tively, illustrating the increasing cosmic shear contribution. The 
right column shows the signal of the same lenses as the middle 
column, but now with misaligned shape models. 

Fig -0 shows the signal of other sets of lens galaxies at 
0.2 < z\ < 0.4 with aligned shape models, namely early types 
at lower stellar mass (9.5 < log 10 M* < 10) in the left column, 
as well as late types with high (10.5 < log 10 M* < 11) and 
low (9.5 < log 10 M* < 10) stellar mass in the middle and right 
columns, respectively. 


4.3.3 Deviations in the isotropic shear profile 

We note that the isotropic AS profile falls below the NFW model 
at very small scales, especially for the lower mass haloes. This 
is caused by the force softening used in the Millennium _Sim- 
ulation and smoothing applied in the ray-tracing dHilbert et alj 
120091: Gillisjh al]l2013h . This is not problematic for our analysis 
given that we only include scales r > 45 kpc//i 73 in the fit. Also, 
the highest stellar mass bins for early types show a marginally 
steeper AS profile than the best-fitting NFW mode l (see Fig.[ 6 ll. 
This could be caused by our use of the iDuffv et all ( 20081) mass- 
concentration relation, which was derived using N-body simula¬ 
tions with a significantly different input cosmology than the one 
used in the Millennium Simulation (<rg = 0.796 versus og = 0.9), 
see e.g. iLudlow et alj d2014l) regarding the cosmology dependence 
of the mass-concentration relation. Since we estimate fh directly 
from the ratio of the measured (/ — / 4 s)AE and AS, the impact 
of these deviations is negligible for our analysis, especially com¬ 
pared to the statistical uncertainties from CFHTLenS. However, as 
a result the model curve for (/ — / 4 s) AS, which is computed from 
the model for AS and the best-fitting /h, is slightly biased high 
compared to the data in Fig.[ 6 ] 


4.3.4 Influence of cosmic shear 

In the two bottom rows of panels in Figs.[ 6 ] and [7] we show the 
anisotropic components /AS and — /45AS, which are decreased 
and increased, respectively, due to the cosmic shear contribution. 
Here the relative effect is stronger towards larger radii, higher red- 
shifts, lower halo masses, and lower /h. In contrast, (/ — / 4 s) AS 
is only weakly affected by the cosmic shear contribution within the 
fitted range (45 kpc//i 7 o < r < r 200 c) as visible in the second row 
of panels, but we note that it is decreased at larger radii compared to 
the case without cosmic shear. This is a result of the non-vanishing 
at these scales (see Sect. ED We stress that our approach using 
the !m 06| formalism and restricting the fit range below r 2 ooc effi¬ 
ciently suppresses the impact of cosmic shear. As visible in Table 
[Dthe remaining net effect is A/h < 0.02 averaged over our red- 
shift range, which is an order of magnitude smaller than current 
observational uncertainties. 


4.3.5 Indications of large-scale structure shape-shear 
correlations 

We note that the anisotropic shear signal shows deviations from the 
simple elliptical NFW model even in the case that no foreground 
cosmic shear is applied to the lens ellipticities, as best visible for 


the early type galaxies in Figs.[ 6 ] and [7] For the early type galax¬ 
ies the measured signal is increased for /AS and decreased for 
—/45 AS compared to the model that assumes an isolated single el¬ 
liptical NFW halo. This trend is the strongest at large radii (for the 
scaling by r shown in the plot), but especially for the lower mass 
haloes it is visible down to small radii (left column of Fig.JD- This 
has the opposite effect compared to the foreground cosmic shear 
which aligns the lens and the source. In fact, both effects actually 
cancel fairly well for the middle column in Fig.[ 6 ] Accordingly, this 
means that the ellipticities of the foreground lenses and the back¬ 
ground sources become anti-aligned because of this effect. This 
is precisely the signature of shape-shear correlations, which are 
one of the intrinsic a l ignment contaminants to cosmic shear (e.g . 


iHirata & Seliakll2004l : Ijoachimi et alj|20nl : iHevmans etal 


201 


and which were studied using the same simulations in i J131 . They 


are caused by alignments of the foreground galaxies with their sur¬ 
rounding large-scale structure, which lenses the background source 
galaxies. Note that the halo ellipticity signal we want to extract 
contributes to shape-shear correlations itself, but an additional con¬ 
tribution comes from the surrounding large-scale structure. We sus¬ 
pect that the additional signal we are detecting here in the context 
of halo shape measurements is caused by this large-scale structure 
contribution. 

Note that this effect appears to be weaker for the 
late type lenses, in agr eement with obse rvatio nal constraints 
dMandelbaum et all 1 201 lh . but also with the |j!3bl intrinsic align¬ 
ments analysis of these simulations. For late type galaxies the 
strongest deviations from the single NFW models are visible for 
the highest stellar mass bin in the middle column of Fig0 Here 
they appear to have the same sign as cosmic shear (alignment of 
lenses and sources), but in absolute terms the effect is weaker and 
the constraints are more noisy than for the early types. 

The observed profiles for /AS, —/ 4 s AS, and (/ — / 4 s) AS 
experience a similar relative suppression when misaligned lenses 
are studied compared to aligned lenses, if foreground cosmic shear 
is not applied (compare the middle and the right column of Fig.[ 6 ]l. 
This is expected given that random misalignment of lens galaxies 
suppresses both the halo shape signature and the additional sig¬ 
nal from shape-shear correlations caused by alignments of the lens 
galaxies with their large-scale environment. In contrast, the appli¬ 
cation of foreground cosmic shear introduces additional additive 
signal which is the same for aligned and misaligned lens models. 

To further test our hypothesis that the extra signal is caused 
by shape-shear correlations from the large-scale environment, we 
further investigate the signal of the aligned early type galaxies at 
0.4 < zi < 0.6 with 10.5 < log 10 M* < 11. For these galaxies the 
plots in the middle column of Fig.[ 6 ] suggest that the contributions 
from cosmic shear and shape-shear correlations cancel approxi¬ 
mately. For these lenses and their surrounding background sources 
we compare the two-point correlation functions 


C± ! (r) = (ei“t7s,t ±ei“*7s,x>(f) 


(25) 


between the (not-lensed) intrinsic lens ellipticities ej nt and the 
background shears (shape-shear signal), versus the corresponding 
cosmic shear correlation functions 


£±° ( r ) = (7i,t7s,t ± 7i, x 7s, x > (r) ■ (26) 

As shown in Fig.[ 8 ] / GI and £ GG are of roughly comparable ampli¬ 
tude (note the opposite sign) over a wide range in scale (r > 0 .' 6 ). 
This is consistent with the interpretation that the influence of cos¬ 
mic shear and the influence of the shape-shear correlations on the 
halo shape measurement cancel approximately. Note however that 


© 2015 RAS, MNRAS POQ.fTHTIl 






























14 Schrabback et al. 




100 1000 
r [kpc/h„] 



100 1000 
r [kpc/h„] 



.1 1 1 

© ® ®_® ®5_® 

*..*.. ******.XI 

.1 . 

CD ® ffi CD ® ® 


.1_ 1 _ j_ 

.1 ' 

100 1000 

L 

1 1 1 ' 1 1 T ^ 

x x x 

.1 - 

-*.**_** * *: 


3 «cTOo'S'0'(d®®®- 0 

■ Ml _ 1 _ 



100 1000 
r [kpc/h„] 


Figure 6. Analysis of the Millennium Simulation showing early type foreground lens galaxies in stellar mass bin 10.5 < log 10 M* < 11: From top to 
bottom we show the isotropic signal AS and the anisotropic signal components (/ — AS, /AS, and AS, respectively. For the anisotropic signal 
components the crosses (circles) correspond to the cases with (without) foreground cosmic shear applied, where the crosses are displayed with an offset 
of half a bin for better readability (note the scaling by r for better readability). The left column corresponds to aligned shape models and lens redshift 
slice 0.2 < z\ < 0.4. The middle column also shows aligned shape models but higher lens redshifts 0.4 < z\ < 0.6. illustrating the effect of an increasing 
cosmic shear contribution. The right column corresponds to the same lens galaxies as the middle column, but now with misaligned shape models, illustrating 
the resulting suppression of the anisotropic signal. The curves indicate the model predictions for the best-fitting /n and best-fitting isotropic model, /n is 
determined from the ratio of (/ — / 4 s) AS and AS for the case without cosmic shear, fitted within the range indicated by the vertical dotted lines. 


mixes contributions from both the halo shape signal (primar¬ 
ily at small r) and the shape-shear correlations from the large-scale 
environment (affecting also larger r). Therefore we do not attempt 
to interpret Fig.[ 8 ]more quantitatively. 

As demonstrated in Sect. l4.3.4l halo shape estimates based on 
(/ — 745 ) AS are insensitive to extra alignment from cosmic shear 
if a limited radial range is included in the fit as done in our analy¬ 
sis. Given that the radial dependence of the measured (/ — / 4 s) AE 
signal is reasonably well described by the model prediction within 
the fit range for most of our lens bins, we suspect that this formal¬ 
ism also provides an approximate correction for the extra alignment 
caused by the surrounding large-scale structure. However, it is clear 
that this model of a single, isolated, elliptical NFW halo is only a 
very crude approximation of reality, and that some deviations due 
to the surrounding large-scale structure are still expected. For ex¬ 
ample, note the significant deviations from the best-fit model also 
within the fit range in the left column of Fig.|T| Nevertheless, as 
we perform a direct, consistent comparison of our CFHTLenS re¬ 
sults to the constraints from the simulation, which ideally includes 
the same physical effects, the relative conclusions should not be 
affected substantially. 


4.3.6 Constraints on fhfor early type galaxies 

We summarise the estimates for the aligned ellipticity ratio /h for 
the different lens samples in Table [3] For early types we generally 
measure larger /h for more massive haloes. For example, the most 
massive low -2 ellipticals yield /h = 0.863 ± 0.011 for the aligned 
lens galaxy models and no cosmic shear. For lower mass haloes 
this value decreases, and we find fh = 0 . 616^0 005 when applying 
cosmic shear and combining both redshift bins and the stellar mass 
range 9.5 < log 10 M* < 11, to best mimic our red galaxy sample 
in CFHTLenS. 

This result is somewhat surprising: Given that the projected 
lens shapes of early types are defined via the projected quadrupole 
tensor of the mass distribution in the simulation, we would naively 
expect to measure /h — 1 in all cases. We suggest three possible 
explanations for this, and we presume that all three may contribute 
in practice at different levels: First, the projected halo shapes are 
approximated via the quadrupole tensor of the mass distribution. 
However, the projected mass distributions will not have exactly el¬ 
liptical shapes, and therefore deviate from the projected light dis¬ 
tribution. Second, the haloes are embedded in a surrounding large- 
scale structure. Components of that large-scale structure that have 
a random orientation with respect to the orientation of the halo will 
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Figure 7. Analysis of the Millennium Simulation showing lenses at 0.2 < z\ < 0.4 with aligned shape models: The left column shows early type galaxies 
with stellar mass 9.5 < log 10 M* < 10. The middle column corresponds to late type galaxies in stellar mass bin 10.5 < log 10 M* < 11. The right column 
also shows the signal of late type galaxies, but with lower stellar masses 9.5 < log 10 M* < 10. For further details see the caption of Fig.jb] 
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Figure 8. Shear-shear and shape-shear correlation functions £ GG and 
—? GI for the aligned early type galaxies in the Millennium Simulation 
at 0.4 < z\ < 0.6 with 10.5 < log 10 M* < 11. For these lenses r 200 c 
equals 0.'98. 


add noise in the form of an on average spherical mass distribu¬ 
tion, reducing fh■ Also the transient sub-haloes, which have been 
removed (see Sect. l4.2l l. contribute to this. Finally, the surrounding 
large-scale structure may have a component that is aligned with re¬ 
spect to the halo orientation, causing the extra signal discussed in 
Sect. l4.3.5l While we expect that our analysis using (/ — / 4 s)AS 
partially separates this component from the signal of the halo, there 
may be some residual effect impacting the fh constraints. A hint for 
this is visible in the left column of Fig.[7] where the model fits the 
measured (/ — / 4 s) AE relatively poorly. We note that the second 
and third effect are likely more pronounced for less massive haloes 
for which the surrounding large-scale structure has a larger relative 
impact, in agreement with the trend observed. 

For the misaligned lens models fh generally decreases as ex¬ 
pected, and we find fh = 0 . 2851 qqq 4 for the combined lens sam¬ 
ple. 

4.3.7 Constraints on fhfor late type galaxies 

For the late type galaxies the recovered values of fh are gener¬ 
ally small, with fh = 0.095 ± 0.005 in the case of aligned models 
and f h = 0.025±q for misaligned models, combining all lens 
samples and applying cosmic shear. This indicates that even if the 
angular momentum vector of each dark matter halo was perfectly 
aligned with the spin axis of its disc galaxy, the halo shape sig¬ 
nal would be largely washed out due to misalignment between the 
mass distribution and the angular momentum vector of the halo. 
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Once galaxy misalignments are included, the net signal is so small 
that it should remain undetected even with much larger surveys than 
CFHTLenS. Despite this generally low signal there is a trend that 
/h increases in the simulated data for more massive haloes. 


4.3.8 Accounting for the differences in the ellipticity dispersions 
between CFHTLenS and the simulation 


We can directly compare our constraints on fh from the Millen¬ 
nium Simulation to those from CFHTLenS if we assume that our 
procedure to assign galaxy ellipticities to the haloes in the simula¬ 
tion is adequate. However, if we compare the ellipticity dispersions 
cr e (including both ellipticity components) of the selected lenses be¬ 
tween CFHTLenS (Table Q} and the simulation (Table [3}. we find 
that the average dispersions are somewhat higher in CFHTLenS, 
with (jCFHTLenS _ q 29 (0 .35) for red (blue) lenses, compared 
to (j| imulation = 0.20 (0.32) for the early (late) types in the Mil¬ 
lennium Simulation. Hence, the simulation appears to deliver on 
average too round lenses, especially for early types. As a first-order 
correction for this discrepancy one could increase the lens elliptic¬ 
ities in the simulation by a factor 


CFHTLenS 

Oe 


^simulation 
u e 


1.46(1.10) 


(27) 


if one assumes that the early (late) type galaxies in the simulation 
resemble the red (blue) galaxies in CFHTLenS. We choose to not 
directly rescale the ellipticities, but to keep the an alysis compara¬ 
ble to earlier studies using these mock catalogues djoachimi et all 
1201 3bl fah. From the definition of fy, (see Eq. I 1 lb it is directly evident 
that such a linear increase of the lens ellipticities e g will reduce fh 
by the same factor given that the shear signal (and hence eh) in the 
simulation is unchanged. Accordingly, when conducting a compar¬ 
ison in Sect.[5]that assumes that the red (blue) CFHTLenS sample 
corresponds to the early (late) type sample in the simulation, we 
will reduce the fh estimate from the simulation by a factor s. 

We also considered alternative approaches to compute scale 
factors s for the approximate matching of the ellipticity dis¬ 
tributions between CFHTLenS and the simulation, such as 
the ratios (|e g |) CFHTLenS /(|e g |) simulation = 1.39 (1.09) or 

«l e g| 2 )/(l e g|» CFHTLenS / «l e g| 2 )/(|e g |» aimulation = 1.52 
(1.11), where the latter is equivalent to matching the (|e g |) with 
an extra weighting by |e g | as employed in Equations and {H>. 
The resulting differences of ±5% (±1%) compared to Eq. ( 127b 
are small, especially compared to the statistical uncertainties from 
CFHTLenS. Also note that we do not need to correct the measured 
lens ellipticity dispersions in CFHTLenS for measurement noise 
given the high signal-to-noise ratio of all lenses in our analysis 7 . 


5 DISCUSSION AND CONCLUSIONS 

Using CFHTLenS observational data we have obtained weak lens- 
ing constraints on /h, the ratio between the aligned projected ellip¬ 
ticity components of galaxy-scale matter haloes and their galaxies. 
In addition, we have used simulated data sets based on the Millen¬ 
nium Simulation to test the methodology and obtain an estimate 
for the expected signal given current alignment models for galaxies 
and their dark matter haloes. 

In this analysis we make use of the formalism introduced by 
lM06l to correct for additional alignments of lenses and sources, 
e.g. due to residual shape systematics or cosmic shear. We demon¬ 
strate that this method removes a simplistic, constant alignment. 


In addition, we demonstrate that a realistic cosmic shear field, 
as present in the Millennium Simulation, is well removed for the 
fit range we employ, leaving residuals A/h < 0.02 for our lens 
redshift range (0.2 < z\ < 0.6), which is an order of magnitude 
smaller than current statistical uncertainties from CFHTLenS. We 
demonstrate that a more simplistic estimator that only considers 
the azimuthal variation in the tangential shear component is highly 
biased (see e.g. the third panel in the rig ht c olumn of Fig .|7]l as 
previo usly pointed out bv lBrainer 5 l2010h and H owell & Brainerd 

boioh . 

From the CFHTLenS data we estimate /h = —0.04 ± 0.25 
for red lenses and /h = 0.69lg’gg for blue lenses combining all 
stellar mass bins with log 10 M, > 10 (log 10 M* > 9.5) for the red 
(blue) lenses. The colour separation is done via the photometric 
type from BPZ, which we use as a proxy for the separation into 
early and late type galaxies. As discussed below, these values are 
broadly consistent with theoretical models. 

In our analysis of the Millennium Simulation we use different 
models for the alignment of galaxies with their dark matter haloes 
for early and late type galaxies. For late type galaxies the analysis 
assumes that the spin vectors of disc galaxies are aligned with the 
angular momentum vectors of their dark matter haloes. In the case 
of perfect alignment we find /h = 0.095 ± 0.005 from the Millen¬ 
nium Simulation, which reduces to fh = 0.025 1° for models 
which assume the misalignment distribution rrom lBettlj20 1 ill . This 
distribution is based on a compilation of results fr om simulatio ns 
that include baryons and galaxy formation physics fo kamoto et alj 
l2005l:ICrain et all2009l:lBett et alJl20ld : lDeason et all201 lh . If we 
rescale these results as discussed in Sect. l4.3.8l to account for the 
differences in the lens ellipticity dispersions in CFHTLenS versus 
the simulation, these values change to /h = 0.086 ± 0.005 for the 
aligned models and fh = 0.0231q'oo 4 for the misaligned models. 
Based on the simulation we do not expect to detect a significant 
signal for late type galaxies in CFHTLenS. Our CFHTLenS con¬ 
straint for blue galaxies is higher than the scaled prediction for the 
aligned models by 1.7cr, and higher by 1.9<r compared to the scaled 
prediction for the misaligned models. 

For early type galaxies the analysis of the simulation assumes 
that the lens ellipticities follow the ellipticity of the projected iner¬ 
tia tensor of the halo mass distribution. Hence, the projected galaxy 
shapes follow approximately the projected dark matter shapes. 
In the case of no misalignment we estimate fh = O. 6 I 6 I 000 ® 
from the simulation when matching the measured halo mass range 
of the CFHTLenS constraints approximately. When applying a 
Gaussian misalignment distribution with an rms scatter of 35° 
(as suggest ed by the distribution of satellites around LRGs in 
SDSS, see lOkumura. Jing & Lil [2009I) . we estimate a value of 
fh = 0.2851q qq 4 . If we rescale these results as discussed in 
Sect. l4.3.8l to account for the differences in the lens ellipticity dis¬ 
persions in CFHTLenS versus the simulation, the values change 
to fh = 0 . 422+0 003 f° r aligned models and fh = 0 . 1951 q ' oo 3 
for the misaligned models. Assuming that the red CFHTLenS 
galaxies directly correspond to the simulated early type galaxies 
would mean that our CFHTLenS constraints are lower than the 
rescaled aligned prediction by 1.8cr, and lower by 0.9cr compared 
to the rescaled misaligned prediction. Hence, they are poorly de¬ 
scribed by the perfectly aligned model but fully consistent with the 
misaligned model. 

It is interesting to compare our con straints to the results of re¬ 
cent hydrodynamical simulations. iTenneti et alj J20I4I) find a mean 
3D misalignment angle at 2 = 0.3 of 25.20° for a halo mass bin 
10 11 ' 5 /i _1 Mq < M < 10 13 /i _1 Mq which most closely matches 
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the mass range of our red galaxies, but this increases to 33.47° 
for lower halo masses 10 lo /i _ 1 M 0 < M < 10 The 
misalignment distribution estimated bv lOkumura. Jing & Li < 2009 1 
and assumed in our simulated analysis (rms s catter of 35°) woul d 
correspond to a similar mean value of ~ 28° iTenneti etall2014h . 
Thus, our assumed misalignment model approximately matches the 
results of this hydrodynamical simulation. However, for the direct 
comparison between observations and simulations there are addi¬ 
tional relevant effects to be considered: for ex ample, baryo ns ap¬ 
pear to make the inner halo more spherical l lBrvan et alj|2013l ). 
In addition, sub-haloes are found to be rounder th an haloes 
iKuhlen, Diemand & Vladaii 120071 : Tenn eti et all 1201 4|) . Accord¬ 
ingly, as our analysis of early type galaxies in the simulation uses 
central haloes only, and given that it does not include baryons, 
we expect that our simulation should somewhat over-predict /h. 
A further complication for the comparison arises from the ob¬ 
servation that numerical simulations of galaxy formation suggest 
that da rk matter halo shapes are m isaligned at d ifferent radii (e.g. 
ISchneider. Frenk & Colell20li IWang et al.ll2014h . 


A number of previous studies have attempted to constrain halo 
ellipticity with weak lensing observationally. Two previous stud¬ 
ies were able to split the lens sample into red and blue galax¬ 
ies as done in our study: lM06l find / h = 0.60 ± 0.38 for red and 
/h = — 1 . 4 + 2 'o f° r bl ue lenses, assuming an elliptical NFW mass 
profile and employing data from the SDSS. Using data from the 
RCS2 and employing the same formalism as I M Obi . I van Uitert et alj 
( 1 201 2 |) find / h = 0 . 20 l}'gj for red lenses and /h = — 2.17+2 03 
for blue lenses when assuming an elliptical NFW mass profile and 
using linear weighting with the lens ellipticity as done in our study. 
While our error-bars appear to be substantially tighter than those of 
these two studies, we note that our analysis of simulated data sug¬ 
gests a sign error in the numerical model prediction computed b y 


pre diction comp uted by 
ytyan Uitert et al. I ( 20121 ) . 


ImOQ for / 45 , which was also employed by[ 

This sign error has likely biased their derived constraints. With the 
corrected sign the predicted signal for (/ — / 4 s)AE is higher at a 
given /h over the entire radial fit range and does not drop as quickly 
towards large r (see Fig.Q}. Accordingly, this correction leads to 
significantly tighter constraints on /h. 

iHoekstra. Yee & Gladdersl l l20041) and Parker et al.1 d2007l) 
use single-band data from RCS and early CFHTLS observa¬ 
tions, respectively, to constrain halo ellipticity without subdivi¬ 
sion into red and blue galaxies. They do not correct for sys- 
tema t ic align men t between lens es and sources as introduced by 
I VI Obi . iHoekstra. Yee & Gladdersl 120041) conduct a maximum like¬ 
lihood analysis assuming elliptical, truncated isothe rmal sphere 
(TIS) models from which they find /h = 0.77 + ° ^. IParker et alj 
d2007l) compute the ratio of the shears measured in quadrants along 
the lens minor and major axes, for which they find a tentative sig¬ 
nal 0.76 db 0.10 when averaged out to 70". There are diff erence s 
in the anisotropic shear fields for NFW and TIS pro files (see MOfj), 
but nonetheless the /h constraint from IHoekstra. Yee & Gladdersl 
d2004l) may appear somewhat high even if their lens selection based 
on magnitude (without colours) would provide a perfect selection 
of early types, which is certainly not the case. The reason for 
this is that their formalism does not account for spurious align¬ 
ment caused e.g. by foreground cosmic shear, which leads to a 
lower measured value of /h. Our analysis of the Millennium Sim¬ 
ulation may provid e a possible exp lan ation for th e high value of 
/h measured by IHoekstra. Yee & Gladdersl d2004l ): As discussed 
in Sect. l4.3.5l we detect an excess signal in /AE for early type 
lens es, which ap proxi mately cor resp onds to the si gnal p robed by 
IHoekstra. Yee & Gladdersl |2004|) and Parker et all 12007 ). This ex¬ 


cess signal has the opposite sign than the signal caused by cosmic 
shear (also visible in —/ 4 s AE). We interpret this signal as the im¬ 
pact of shape-shear intrinsic alignments, whi ch are a major con¬ 
taminant to cosmic shear measurements (e.g . lHirata & Seliakl2004l : 
Ijoachimi et alj l20lTI : iHevmans et all |2013h . They are caused by 
an alignment of foreground galaxies with their surrounding large- 
scal e structure (as e.g. detected in the distribution of red galaxies , 
see iMandelbaum et al1l2006bl : iLi et alj 120131 : IZhang et al.l |2013h . 
which also lenses the background sources. While the halo elliptic- 
ity signal co ntribu tes to small-scale shape-shear correlations itself 
jBridle & Abdallall2007l) . there appears to be an additional compo¬ 
nent generated by the large-scale structure the halo is embedded in. 
In other words, the assumed model of an isolated elliptical NFW 
halo is too simplistic. The net effect of this signal is a net anti¬ 
alignment (orthogonal alignment) of the ellipticities of foreground 
lenses and background sources. We expect that the formalism in¬ 
troduced by iMQq and employed by us also provides a partial cor¬ 
rection for this contaminant, but conclude that the most robust halo 
shape results can be obtained from the direct, relative comparison 
between observations and simulations that include this effect, as 
done in our study. 

We note that our analysis of blue galaxies in CFHTLenS 
shows indications for the influence of cosmic shear on the (un¬ 
corrected) anisotropic shear components /AE (slightly negative 
signal in the third row of panels in Fig.0 and — /45AE (slightly 
positive signal in the fourth row of panels in Fig. 0. In contrast, 
we do not observe this trend for red galaxies (see Fig.[3]i. A likely 
explanation for this may be that the additional shape-shear signal 
discussed above appears to roughly cancel the cosmic shear con¬ 
tribution for our red galax y sample. In cont rast, the magni tude- 
selected le ns samp les from lHoekstra. Yee & Gladdersl J20041) and 
IParker et al J ( 12007 3 will likely be dominated by galaxies at lower 
redshift. Here, the shape-shear contribution will dominate given the 
smaller cosmic shear signal. 

As discussed in Sect. l4.3.6l the influence of neighbouring 
large-scale structure is likely one of the main reasons why the anal¬ 
ysis of the aligned early type galaxies in the Millennium Simulation 
yields somewhat smaller /h than the naive expectation of /h — 1 , 
in particular for the lower mass haloes. 

The early type lenses with the highest halo mass show the 
largest /h in our simulation (see Table[3]i. In ad dition, more massive 
haloe s are less spherical in simulations (e.g. iBailin & Stein met/1 
1 20051 : iDespali, Giocoli & Tormenl 120141) . Furthermore, galaxies 
with more massive haloes are expected to be less misaligned 
iTenneti et alj|2014l ). From all three effects we expect that future 
studies with larger samples might have the best prospects for de¬ 
tecting halo ellipticity with weak l ens ing for very massive early 
type galaxies. Very recently, IClampitt & Jainl J2015I ) reported the 
detection of a significant halo shape signal for exactly such lurni- 
n ous red galaxies from SDSS, employing a new estimator from 
lAdhikari. Chue & Dalall i2015l) . In the future, it will be interest¬ 
ing to test this estimator on large mock data-sets, such as the data 
provided by the Millennium Simulation. 

Weak lensing studies constraining halo ellipticity are not 
only interesting for a better understanding of the link between 
galaxies and their surrounding matter haloes, but have also 
been discussed as a possible test for theories of modi fied grav¬ 
ity su ch as MOND (Modified Newtonia n Dynamics. iMilgroml 
119831) . TeVeS (Scalar-Tensor-Vector theorv. [Bekensteinl2004 ), and 
MOG/STV G (Modified Gravity and Sca l ar-Te nsor-Vector Grav¬ 
ity theory, iMoffaj 120061 : iMoffat & Tothl I2009I). for which lens¬ 
ing prescriptions have been developed iMortlock & Turneill200 ll : 
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iBekensteinl |2004 IChiu. Ko & fiaiil 120061 : iMoffat & Tothl l2009h . 
The equival ent signa l s to halo sh apes are discussed for isolated 
galaxies in iMilgroml J200lh and Isellwood & Kosowskvl l l2002h . 
predicting an isotropic shear signal towards large radii. Accord¬ 
ingly, a significant detection of halo ellipticity from weak lens- 
ing, which would be expected within ACDM for the analysis of 
early type galaxies in upcoming experiments, could be interpreted 
as evidence against such theories of modified gravity. However, our 
analysis of the simulated data has shown that already within ACDM 
neighbouring structures have a considerable influence, which needs 
to be taken into account. Thus, better model predictions need to 
be developed also for theories of modified gravity, that include 
both large-scale structure contributions and the influence of baryon 
physics. 

The prospects for near-future improvements on halo shape 
constraints from weak lensing are relatively good: On the one 
hand, weak lensing surveys are improving rapidly in size. Surveys 
w hich are already und erway include th e Dark Energy Survey (DES, 
iThe Dark Energy Survey Collaborationl2005)), the Hyp er S uprime - 

Cam Subaru Strategic Program (H SC-SSP.lM ivazaki et ~ai1l2012h . 

the Kilo Degree Survey (KiDS, Ide Jong et alJ l2013h . as well 
as the Panoramic Survey Teles cope and Rapid Response Sys¬ 
tem (PanSTARRS, Kaiser et~ah| 20ic|). and future pro grammes 
such as LSST JLSST Science Collaboration et alj|2009h and Eu¬ 
clid l lLaureiis et al J 2011) will tighten parameter constraints even 
further. On the other hand, the use o f new estimators may 
also provide additional insight: Simon. Schneider & Ktibleii 120121) 


and [Adhikari. Chue & Dalall d2015l) suggest that the position an¬ 
gle dependence of the galaxy-shear-shear correlation function, 
which is one of the obse rvables of galaxy-galaxy-galaxy lensing 
dSchn eider & Wattsl l200^) . may provide constraints on halo ellip¬ 
ticity independent of the orientation of the lens galaxy ellipticity. 
This measurement would be unaffected by misalignment, but is 
expected to have lower signal-to-noise. In addition, the measure¬ 
ment of higher-order lensing (“flexion”) may als o provide addi¬ 
tional sensitivity for constrai ning halo ellipticity dEr & Schneideil 
|201 ll : Er & Bar telmamill 20131) . 
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APPENDIX A: TESTS WITH THE “CLONE” IMAGE 
SIMULATION 

iMiller et all<l2013h and lHevmans et alJd2012h present detailed tests 
of the lensfit shape estimation algorithm for cosmic shear measure¬ 
ments. As the required level of systematics control is less demand¬ 
ing for galaxy-galaxy lensing studies than for cosmic shear, we do 
not present general shape measurement tests again in the current 
work. However, what has not been tested before in detail, is the 
recovery of galaxy shapes in the presence of a nearby bright lens 
galaxy, whose light might affect the shape measurement process 
for sources at small angular separations (e.g. smaller than the 9" 
postage stamp sizej3 T° test if this could be a concern for our 
study, we investiga te len sfit sh ape m easurements of the “Clone” 
image simulations IMiller et al]|2013l) . These simulations contain 
simulated galaxy images with the positions and magnitudes from 
CFHTLenS and cosmological shear from an N-body simulation 
dHarnois-Deraps. Vafaei & Van Waerbeke! 120121) as input, but no 
galaxy-galaxy lensing signal. Hence, any detected galaxy-galaxy 
lensing signal would indicate a spurious effect introduced by the 
measurement process. 

The result of this test is shown in Fig. lAll Here we split 
the “lenses” into magnitude bins and have subtracted the intrin¬ 
sic ellipticity and cosmological shear from the measured elliptici- 
ties to maximise the S/N of the test and achieve statistical error- 
bars which are significantly smaller than for the actual survetF^l 
Given the lack of a clear signal at small scales, and the significantly 

14 For example, a circular source might appear slightly elliptical due to 
contamination by light from the nearby lens, where the net ellipticity would 
point towards the lens, thus generating spurious signal with AS < 0. This 
effect would likely be stronger along the direction of the lens major axis, 
and could thus al so generate / AS < 0. 

15 In contrast tolMille r et all 120131) . we do not make use of rotated galaxy 
pairs in the simulations, which are often used to partially cancel shape noise 
from the intrinsic source ellipticities. This combination is not useful for our 
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Figure Al. Test for a potential systematic contamination of the small-scale halo shape signal originating from lensfA shape measurements in the presence of a 
nearby bright galaxy (lens): We employ lensfit shape measurements of the “Clone” image simulation, which uses actual galaxy positions and magnitudes from 
CFHTLenS as input but no galaxy-galaxy lensing signal. To optimise the sensitivity of the test, all bright galaxies in the indicated i- band magnitude intervals 
with lensfit shape estimates in the CFHTLenS data were used as lens sample (left: 18 < i < 19, middle: 19 < i < 20, right: 20 < i < 21), irrespectively of 
their assigned CFHTLenS redshift (we assume our central lens redshift z = 0.4 for the computation of angular diameter distances). The top (bottom) panels 
show the isotropic (anisotropic) shear signal, where both have been scaled by r for better readability of the plot. 


smaller uncertainties compared to the CFHTLenS constraints, we 
conclude that the light of nearby bright lenses does not appear to 
introduce significant spurious signal for our halo shape analysis. 


test, as also the nearby “lenses” will be rotated, therefore leading to different 
lens-light contributions in the rotated source galaxy pairs. 
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